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Abstract 



The scale dependence of an effective average action for mesons and quarks is 
described by a nonperturbative flow equation. The running couplings lead to spon- 
^ ■ taneous chiral symmetry breaking. We argue that for strong Yukawa coupling be- 

tween quarks and mesons the low momentum physics is essentially determined by 
infrared fixed points. This allows us to establish relations between various param- 
eters related to the meson potential. The results for f n and (^ip) are n °t very 
sensitive to the poorly known details of the quark-meson effective action at scales 
where the mesonic bound states form. For realistic constituent quark masses we 
find f n around 100 MeV. 
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1 Introduction 



Quantum chromodynamics as the theory of strong interactions and its symmetries are well 
tested both for high and low momenta. For momenta q 2 ^ (2 GeV) 2 asymptotic freedom 
[1] permits the use of perturbation theory for a quark-gluon description with small gauge 
coupling g s . The long distance behavior for q 2 ^ (300 MeV) 2 can partially be described 
by chiral perturbation theory [2, 3]. Here the picture is based on a nonlinear or linear 
cr-model [4] for the pseudo-scalar mesons. The latter can also be extended to describe in 
addition scalar, vector and pseudo-vector mesons. Such effective models incorporate the 
chiral symmetries of QCD and use several free couplings to parameterize the unknown 
strong interaction dynamics. The parameters are determined phenomenologically [3, 5], 
but on a more fundamental level the question arises how they can be related to the 

2 

parameters of short distance QCD, i.e. the strong fine structure constant a s = ^ and 
the current quark masses. For example, one may ask how the most prominent quantity 
of the mesonic picture, namely the pion decay constant /„- which measures the strength 
of spontaneous chiral symmetry breaking, can be computed from a s or vice versa. 

Important progress in this question has been achieved by numerical simulations in 
lattice gauge theories [6]. Serious difficulties in this approach remain, however, related 
to the treatment of dynamical quarks and chiral symmetry. There is also a vast amount 
of literature on various analytical attempts to attack this problem. Examples can be 
found in [7]. In this paper we employ a new analytical method based on nonperturbative 
flow equations for scale dependent effective couplings. These couplings parameterize the 
effective average action T k [8] which is a type of coarse grained free energy. It includes 
the effects of all quantum fluctuations with momenta larger than an infrared cutoff ~ k. 
In the limit where the average scale k tends to zero Tk-*o becomes therefore the usual 
effective action, i.e. the generating functional of 1PI Green functions [9]. The scale 
dependence of T k can be described by an exact nonperturbative evolution equation [10, 9] 

|r tM = lT r {(r« 2 »M + « t )-'^} (i.i) 

where t = ln(/c/A) with A some suitable high momentum scale. The trace represents here 
a momentum integration as well as a summation over internal indices and we note the 

(2) 

appearance on the right hand side of the exact inverse propagator T k as given by the 
second functional variation of with respect to the field variables ip. The function R k (q) 
parameterizes the detailed form of the infrared cutoff or the averaging procedure. With 
the choice 1 

Z ua 2 e- q2 / k2 

*M> = t-V-^ < h2 > 

we observe that the momentum integration in eq. (1.1) is both infrared and ultraviolet 
finite. For fluctuations with small momenta q 2 <C k 2 the infrared cutoff R k ~ Z^ jk k 2 
acts like an additional mass term in the propagator, whereas for q 2 ^> k 2 it is ineffective. 

1 Z Vy k is an appropriate wave function renormalization constant which will be specified later. 
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The only difference between the flow equation (1.1) and the /c-derivative of a one-loop 
expression with infrared cutoff R k concerns the appearance of instead of the second 
functional derivative of the classical action. This turns eq. (1.1) into an exact equation, 
but also transmutes it into a complicated functional differential equation which can only 
be solved approximately by truncating the most general form of IV As it should be, 
eq. (1.1) can be shown [11] to be equivalent to earlier forms of the exact renormalization 
group equation [12]. It may be interpreted as a differential form of the Schwinger-Dyson 
equations [13]. The difficult part is, however, not so much the establishment of an exact 
flow equation but rather the finding of a suitable nonperturbative truncation scheme 
which allows to solve the differential equation. Then the flow equation can be integrated 
from some short distance scale A, where I\ can be taken as the classical action, to k — > 
thus solving the model approximately. 

Within the formalism of exact flow equations for the average action it is possible to 
change the relevant degrees of freedom [14]. The idea is now to start from the exact flow 
equations for quarks and gluons for k > k v , and to use a similar exact flow equation 
for quarks and mesons for k < k v . The transition at the scale k v (600 — 700 MeV) 
between the two pictures can be encoded into an exact identity [14] which replaces multi- 
quark interactions in the quark-gluon picture by mesonic interactions in the quark-meson 
picture. We emphasize that in the quark-meson description the quarks remain important 
degrees of freedom as long as k is larger than a typical constituent quark mass m q ~ 
300 MeV. We have therefore to deal with an effective quark-meson model, where the 
mesons are described by a linear cr-model with Yukawa coupling h to the quarks. A 
first attempt to describe the transition to a quark-meson model within a QCD-inspired 
model with four-quark interactions [14] has been very encouraging. Spontaneous chiral 
symmetry breaking was observed for low k, with a chiral condensate of the right order 
of magnitude. In the following, the formalism has been generalized to QCD [15], with 
a method where the gluonic fluctuations with q 2 > k 2 are integrated out subsequently 
as k is lowered. So far, the treatment of the quark-meson model for scales k < k^ has 
been very rough, however, since only quark fluctuations were included in ref. [14]. It is 
the purpose of this paper to present a systematic study of the scale dependence of the 
quark-meson effective action, including both quark and meson fluctuations. 

Our main tool are nonperturbative flow equations which describe the change of shape 
of the effective meson potential and the running of the Yukawa coupling h(k). In the 
perturbative limit these equations reproduce the running of the couplings [16, 17] in a 
Ul{N) x Ur(N) model, which has been investigated in the context of dynamical top 
quark condensation [18]. In our context N stands for the number of quark flavors. The 
most important nonperturbative ingredient in the flow equations will turn out to be the 
appearance of effective mass threshold functions which account for the decoupling of 
modes with mass larger than k. The solution of our approximate flow equations allows 
us to express typical low momentum quantities like f n in terms of the "initial values" for 
the quark-meson model at the scale k v , as for example the meson mass term Tn 2 (k ip ) or 
the wave function renormalization constant Z^ik^). 

Not too surprisingly, the effective quark-meson Yukawa coupling h will turn out to be 
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rather strong. This can easily be seen by noting that for k — > this coupling is related 
to the ratio of a constituent quark mass m q to f n = 93 MeV, namely 



h(k = 0) 



2m q 



~ 6.5 . 



(1.3) 



For larger k the Yukawa coupling must be even stronger, with a typical nonperturbative 
initial value h 2 {k v ) /167r 2 <; 1. The presence of a strong coupling has important conse- 
quences for the predictive power of the quark-meson model. Generically, the system of 
flow equations exhibits (partial) infrared fixed points in the absence of a mass scale. Due 
to the large Yukawa interaction the couplings are driven very fast towards these fixed 
points and the system "looses its memory" on the detailed form of the initial values at 
k v . Despite the fact that the running is finally stopped by the formation of the chiral 
condensate this infrared stability implies that the low momentum quantities essentially 
depend only on one "relevant" parameter at the scale k v , i.e. the ratio m 2 (/c v )/A; 2 . Using 
the value (1.3) for h(0) and h(k v ) between 12 and 100 we find in a simplified model 



in good agreement with the observed value f n = 93 MeV. An estimate of the error as well 
as a more complete treatment including more accurately the effects of the chiral anomaly 
and the strange quark mass is postponed to future work. 

Besides the exciting prospect of computing f n and other parameters of the low mo- 
mentum meson interactions from QCD our approach seems also capable to deal with 
other issues. Once the parameters at the transition scale k v are fixed either by a QCD- 
computation or by fitting low momentum observational data, it is straightforward to study 
the quark-meson system at nonvanishing temperature. With methods described in ref. 
[19] the temperature dependence of f n or (i/jip) can be investigated. For T ^ k^/2^ a 
study within the effective quark-meson model with T-independent initial values at k v 
should be sufficient, whereas for larger temperatures the T-dependence of initial param- 
eters like rn 2 {k v ) starts to become an important effect. One may therefore hope to gain 
new insight into the nature of the chiral phase transition in QCD [20]. Another inter- 
esting issue concerns the use of quark-meson models to describe hadronization in high 
energy scattering experiments involving quarks or gluons [21]. Here our approach may 
help to compute the phenomenological parameters used in those models. Finally, the 
scalar-fermion models have been extensively studied in the large-iV c limit [17, 18], and 
our nonperturbative flow equations may help to access smaller values of iV c . 

Our paper is structured as follows: in section 2 we give a brief phenomenological 
introduction to the chiral quark-meson model with N flavors. The scale dependence of 
the effective meson potential is then described in section 3 and the scalar wave function 
renormalization can be found in section 4. In section 5 we derive the /9-function for the 
running Yukawa coupling between quarks and mesons as well as the quark wave function 
renormalization. Section 6 is devoted to a short discussion of the chiral anomaly and the 
presentation of two simplified models with two quark flavors. The first model is based 



U - (83 - 100) MeV 



(1.4) 
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on the symmetry C/l(2) x Ur{2) and neglects the effects of the chiral anomaly, whereas 
the second one based on 0(4) neglects all scalars whose masses obtain contributions from 
the chiral anomaly. In section 7 we discuss in detail the infrared stability properties 
for models with strong Yukawa couplings and the consequences for the "prediction" of 
/„-. Section 8 finally contains our quantitative estimates for f n and the chiral condensate 
(xjjxjj^. Conclusions are drawn in section 9. 



2 The chiral quark— meson model 



We describe the low-energy degrees of freedom of QCD by an effective action for quarks 
and mesons. We concentrate in this paper on pseudo-scalar and scalar mesons if which 
transform in the (N, N) representation of the flavor symmetry group SUl(N) x SUr(N) 
for N flavors. We consider the chiral limit where the current quark masses are neglected. 
In its simplest form the effective action T k for quarks and mesons contains kinetic terms, 
a potential for the scalar fields and a Yukawa coupling between quarks and mesons: 



r fc = J 'd 4 x{z ip (k)d lx < P * ab d^ ab + u k (<p,<p*) 



(2.1) 



Our Euclidean conventions (h(k) is real) are specified in appendix A. The scalar potential 
is assumed to be a function of the invariants 



P 



tr (</?V) 

N ( t V 

det ip + det ^ 



N — 1 



(2.2) 



where we neglect the dependence on additional higher order invariants present for iV > 3 
(cf. appendix B). 

Spontaneous chiral symmetry breaking with a residual vector-like SU(N) flavor sym- 
metry occurs if the potential has a minimum for a ^ 



V 



\ 



00 / 



Po = N \a \ 



(2.3) 



In this case we consider a quartic approximation for the potential 

1- N — 1— 1 

U k = -JI 2 (k)p + -M^p 2 + -^\ 2 (k)r 2 - -V(k)£ 



(2.4) 
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where p 2 {k) is related to the /c-dependent minimum value po(k) by 



N-2 
2 



^ k)= j M) -m{mj . ,2.5) 

Without loss of generality we will restrict ourselves to positive V. Up to an irrelevant 
constant we can also write 

N-2 

^^(p-po^V^A^ 2 p. (2.6) 

On the other hand, at short distance scales spontaneous chiral symmetry breaking is not 
yet visible and £4 is in the symmetric regime (a = 0), where we use the parameterization 

U k = m\k)p + h^p 2 + ^^A 2 (A;)r 2 - l -v{k)i . (2.7) 

With these approximations our model can be described in terms of the renormalized 
couplings 

h(k) = Z-V\k)Z^{k)h{k) 
A 1>2 (*) = Z~ 2 (k)\ 2 {k) 

and either the mass term 

m 2 (k) = Z^ 1 {k)m 2 (k) (2.9) 
or the location of the potential minimum 

p R (k) = Z 9 {k)p Q {k) . (2.10) 

The dimension of V depends on iV and the renormalized coupling is 

_ AT 

u R (k) = Z v 2 (k)u(k) . (2.11) 

For 17 = the model has an additional axial Ua(X) symmetry which, however, is broken 
in QCD through the axial anomaly. 

The quark-meson model is supposed to be obtained as an effective model at some 
scale ky, say ~ 600 MeV. It should be derivable from QCD by integrating out the 
gluonic degrees of freedom and converting nonlocal four-quark interactions into an ef- 
fective quark-meson theory by the change of variables described in [14]. This gives a 
direct relation between ip and a suitably defined [14] composite quark bilinear operator 

O a b = ^Va) according to 

b _ 2m 2 (k ip )Z il (k ip ) b 

The aim of this paper is to follow the evolution of from the "initial value" at k = k^ 
to k — 0. The effective action T = Tk=o then describes the 1PI Green functions for the 
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collective meson fields or quark bilinears. In particular, the chiral condensate is related 
to the vacuum expectation value of ip corresponding to the minimum of the effective 
potential U = U k =o through 2 

with 3 

ffo= (*fco)p ( , 14) 

For o"o different from zero the chiral symmetry is spontaneously broken and the spec- 
trum of scalars contains N 2 — 1 Goldstone bosons corresponding to the pions. Their 
interactions are described by the nonlinear cr-model. Neglecting the explicit SUv(N) 
breaking through quark masses the pion decay constant f n is given in our conventions by 

U = 2a R (2.15) 

with renormalized expectation value 

a R = z'j 2 (k = o)<r„ = j 1/2 . (2 , 6) 

In our normalization the experimental value reads = 93 MeV or 

a R = 46.5 MeV . (2.17) 

Another interesting quantity in our picture is the renormalized quark mass (h = h(0)) 

m q = h(jR (2-18) 

This corresponds to a constituent mass generated by chiral symmetry breaking. (We 
remind that we consider the approximation of vanishing current quark masses here.) A 
typical value should be around 300 MeV and the renormalized Yukawa coupling therefore 
be relatively large, h ~ 6.5. 

The scalar spectrum is discussed in detail in appendix B. For a^^O the meson sector 
contains besides the iV 2 — 1 massless Goldstone bosons the cr-field (radial mode). With 
Ai = Ai(0) and vr = v>r(0) its mass is given by 

ml = 2iVA 1( r 2 - l -v R {N - 2)<~ 2 . (2.19) 

One meson acquires a mass through the chiral anomaly. For the realistic case of = 3 
this can be identified with the r]' meson whereas the 77 meson remains massless in the 

2 We note that (2.13) defines the quark condensate at the scale k v which may be different from the 
scale of usual chiral perturbation theory estimates. 

3 Neglecting fermion masses the phase of do is arbitrary and we employ here a real and positive 00. 
Correspondingly, (ip^) stands only for the magnitude of the chiral condensate. 
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chiral limit as one of the Goldstone bosons. For N = 2 we are left with the three pions 
as massless degrees of freedom whereas the if-mesons and the 77-meson are absent from 
the spectrum. We will also for N = 2 associate the anomalously massive meson with the 
rj' meson. In our model its mass is given by 

™ 2 n , = y^£" 2 • (2.20) 

The remaining N 2 — 1 massive scalar fields in the adjoint representation of the diagonal 
flavor symmetry group SU(N) have mass (for A2 = ^(0)) 

m 2 a = N\ 2 a 2 R + v R a%- 2 . (2.21) 

The neutral component can be associated with the ao meson with mass 983 MeV. For 
realistic meson masses (m^ = 958 MeV) the couplings would be 

v R ~ (958 MeV) 2 for N = 2 

v R ~ 13158 MeV for N = 3 

A 2 ~ 11 for N = 2 

(2 23) 

A 2 ~ 55 for N = 3. 

One should, however, notice that the values for u R and A 2 are rather sensitive to the 
precise association of m a or m v / with known particle masses and should therefore only be 
taken as a rough estimate. In particular, for N = 3 the effects of a nonzero strange quark 
mass have to be incorporated for a more realistic estimate. 

3 Scale dependence of the effective meson potential 

The meson degrees of freedom can be introduced [14] at some short distance scale k v by 
inserting the identity 

1 = const J Va A Va H exp j-1 [(a\ - K\<G - Oty]G) CT 1 (a A - GK A - GO 
+ (4 - KjjG - 0^]G) CT 1 (a H - GK H - GO®[il>])] } 

(3.1) 

into the functional integral for the effective average action for quarks. Here K a ,h are 
sources for the collective fields and correspond to the antihermitian and hermitian parts 4 
of (p. They are associated to the fermion bilinear operators 0[ip], 0^[ip] whose Fourier 
components read 



— ^g(-p,p + q)ip a (p)MP + q) 
J p _ a (3-2) 

o {b) \{q)=- j j^^g(-p,p + q)4> a (p)iMp + q) ■ 



4 The fields a ah associated to Kah by a Legendre transformation obey a a — —^(cp — ^), &h = 
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The wave function renormalization g(—p,p + q) and the propagator G(q) are chosen such 
that the four-quark interaction contained in (3.1) cancels the dominant part of the QCD- 
induced nonlocal four-quark interaction in the effective average action formulated only 
for quarks. As a result, the introduction of collective fields by (3.1) replaces the dominant 
part of the four-quark interaction by terms quadratic and linear in the meson field. The 
resulting effective quark-meson interactions are more general than those of the model 
described in the last section. The momentum dependence of the kinetic terms and the 
Yukawa couplings can be described by an extended truncation of the effective average 
action which, for general space-time dimensions d, is given by 

r fc = Jd d xu k ( v ,^) 

+ / (^rf { Z ft^)? 2tr (^(^(?)) + ^,fc(#(?)7"?^(?) (3.3) 

+ / -^^k{-q,q-p)${q) (j-^-<p(p) ~ ^^V(-p)) • 

At the scale k v the average potential is then purely quadratic 

U kip =m 2 tr(^V) (3.4) 

and the inverse scalar propagator is related to G(q) in eq. (3.1) by 

G~\q) = 2m 2 + 2^(g)g 2 

— (3.5J 
Z v (q) = Z vK (q) . 

The initial value of the Yukawa coupling corresponds to the "quark wave function in the 
meson" in eq. (3.1), i.e. 

h kv (-q,q-p) = g(-q,q~p) (3.6) 

which can be normalized with ^(0,0) = <?(0, 0) = 1. The propagator G and the wave 
function g(—q, q—p) should be optimized for a most complete elimination of terms quartic 
in the quark fields. Neglecting the remaining if) A terms and terms of higher order in ip 
(e.g, ip 6 ) one arrives at the initial value for 1^ . In the present paper we often do not want 
to keep the complete momentum dependence of Z^ tk , Z Vtk and h k . Useful definitions of 
the initial values of the parameters of the model in section 2 are then 



q 2 =0 



Z^ky) = Z^ kip (q) 
Z*{K) = ^{G-\q)-G-\Q)) 
m\k,) = \g-\Q) 



(3.7) 



2 



KK) = h k v (-q,q) 



= i . 

q 2 =0 
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Although the results of [14] should only be considered as rough estimates it seems 
convenient to use them as a guide for the choice of initial values of the various couplings. 
The values found for the transition scale k v and the scalar mass at this scale are [14] 
k v = 630 MeV, m(k v ) = 120 MeV. The g 2 -dependence of G was not computed very 
reliably in ref [14]. Large- A^ estimates use a rather weak g 2 -dependence [17]. As a typical 
guess we consider here, somewhat arbitrarily, that G~ 1 (q 2 = /c 2 ) exceeds G~ 1 (q 2 = 0) by 

15%. This leads to Z v (k v ) = 0.15^^ ~ With Z^[k v ) = 1, h(k v ) = 1 this 

corresponds to a large renormalized Yukawa coupling of h 2 (k v ) = 180. We will see later 
(section 7) that for strong initial Yukawa couplings the decisive parameter is the ratio 

to = Z}(K)^^. (3.8) 

<f 

The values of [14] correspond to eo = 0.036. Both, Z^ and m, may be somewhat lower 
than the values from [14] and we will often use typical values mZ^(k v ) = 89 MeV(63 MeV) 
for which e = m 2 Z 2 (k v )/k 2 ~ 0.02(0.01) and h 2 (k v ) ~ 330(660) if the same assumption 
on the momentum dependence of G~ l (q) is made as above. The dependence of our results 
on the choice of initial values will be discussed in detail in section 7. 

We note that the use of the identity (3.1) does not lead to anomalous ^7^(1) violating 
meson interactions and (3.3) conserves the axial ^7a(1) symmetry. Consequently the solu- 
tion of the flow equations for the /c-dependence of the average potential also conserves this 
symmetry. The formalism has therefore to be extended to incorporate anomalous fermion 
interactions of the type Det (^jj a ipbj into the mesonic picture. This issue may be addressed 
for the time being by introducing a term —\vt; into C/ fc of (3.4) as a phenomenologically 
determined coupling. We leave this for future work and concentrate here on the Ua(1) 
conserving case V = and later (section 6) on the opposite extreme u — > oo for N = 2. 
We also observe that a generalization of the formalism of [14] may lead to nonvanishing 
meson self-interactions Ai, A 2 at the scale k v as predicted by large- N c results [17]. 

At the scale k^ the effective potential (3.4) has its minimum at the origin. As a result 
of quantum fluctuations with momenta q 2 < k 2 one expects that the potential changes its 
shape and ends up at k = with a minimum for p > 0, resulting in a spontaneous breaking 
of chiral symmetry. The aim of this paper is to derive flow equations for the /c-dependence 
of Zijj, Z<p, Uk and h and to compute the observable quantities at k — described in the 
last section from the initial values (3.7). Solving the flow equations numerically we find 
that chiral symmetry breaking indeed occurs as demonstrated in figure 1. 

We begin with the evolution equation for the effective average potential Uk- Except 
for the /c-dependence of V the evolution equation for the potential can be obtained (c.f. 
appendix A) by studying a constant scalar field configuration which is real and diagonal 

Vab = VaSab = m a 5 ab (3.9) 

We evaluate the exact evolution equation (1.1) for this configuration and insert the trun- 
cation (2.1) into the right hand side. The flow equation has a fermionic and a bosonic 
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Figure 1: Evolution of the renormalized mass m in the symmetric regime (dashed line) and the vacuum 
expectation value an of the scalar field in the SSB regime (solid line) as functions of k for the Ul(2)x Ur(2) 
model. Initial values are Xi(k lp ) — \2{k v ) = for k v = 630 MeV with h 2 {k v ) — 300 and e = 0.01. 



contribution. The bosonic part follows by neglecting for a moment the quarks [10, 9]: 



d_ 

dt 



U, 



1 r d d q d 



2 J (27i) d dt 



i 



+ 



i 



+ E 



a L Z^P + Mi, Z v P + M? a 
1 



+ 



Z V P + (M+J* z v p + (M^ b y 



+ 



Z v P + (M+ b y Z v P+(M Fab )\ 
We observe the appearance of the (massless) inverse average propagator 



(3.10) 



P{q) =q 2 + Z- x R k {q) 



1 - e ^ 



(3.H) 



which incorporates the infrared cutoff function R k (1.2). The dependence of the various 
mass eigenvalues on tp a can be found in appendix B. We restrict the discussion here to 
the approximation V = which corresponds to neglecting the mass difference between 
the pseudo-scalar pion triplet and the rj singlet. The mass eigenvalues on the right hand 
side of (3.10) are then given by (B.10) — (B.15). In order to express the eigenvalues m 2 a in 
terms of p and r 2 we consider a particular configuration ip where N — 1 eigenvalues are 
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equal to m\ such that 

p = (N — l)m\ + m 2 N 
tr0 = N (^i - m N j , (f) = <pi<p-—p 

or 

1 (3.13) 

This defines the right hand side of the evolution equation (3.10) as a function of p and 
r 2 . In the symmetric regime the evolution of the couplings Ai, A 2 and the mass term m 2 
can now be extracted by suitable differentiation of eq. (3.10) with respect to p and r 2 , 
evaluated for p = r 2 = 0. One finds for the bosonic contributions 

d_ 2 1 r d d q dR k 2{N 2 + l)Ai + {N 2 - 1)A 2 





N 




- — -2 




= m\ 








= ™?n 



-m 



1 r d d g dR k 2(N' + l)X 1 + (N'-l)X 2 
dt"" ~ 2 J (27i) d dt {Z v P + m 2 f 1 ' 

d_j f d d q dR k 2(N 2 + A)\\ + 2(N 2 -l)\{\ 2 + (N 2 -l)\l 

dt 1 7 (2vr) d dt (Z^P + m 2 ) 3 1 j 

0_ _ /• dR k 12\ 1 \ 2 + 2(N 2 -3)X 

Ot 2 y (2ir) d dt (Z^P + m 2 ) 3 ' 1 J 

If in the course of the evolution towards smaller values of k the mass term m 2 becomes 
negative we should switch to the couplings appropriate to the regime with spontaneous 
symmetry breaking (SSB regime). There we define 

Ai = U£(po,T 2 = 0) 

T 4 du *( m (3 ' 17) 

= iv3TW Po ' r2 = 0) 

where p corresponds to the fc-dependent minimum of the potential. We use that 

U'M = (3.18) 

is valid for all k and therefore obtain 

The evolution equations for p , Ai and A 2 follow directly from the definitions (3.17), (3.19). 
For V = they read 



a i r d d q dRk { n 2 3 _ + t 

9r° 27 {2n) d dt I (Z V P) 2 ^ (^P + 2A lPo ) 2 (^P + A 2 p ) 9 ' " 
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d_ 
di 



Ai 



d d q dR, 



9 \ 
d~t 2 




9X1 



z^p + 2\ lPo y 



+ 



(N 2 -l) (A! + A 2 ) 2 > 
(^P + A 2Po ) 3 



(3.21) 



+ 



+ 



(27r) d dt 
9(N 2 - 4) 



(V) J 



A 



+ 



iV 2 A 2 



[Z V P + A 2Po ) 3 4p, 



3A 2 (|A 2 + A!) 
M5A2 - Ai) 



:^p + a 2Po ) 2 (z v py 



1 



(3.22) 



;^P + 2A lPo ) 2 (^P + A 2Po ) 5 



It is straightforward to check that in the limits m 2 — > 0, po — y the flow equations for Ai, 
A 2 coincide in the symmetric and SSB regime. We also note that the evolution equation 
for A 2 depends on the precise definition of this coupling. This issue is shortly addressed 
in appendix C where we also give an alternative formulation of eq. (3.22). 

Next we turn to the fermionic contribution to the evolution equation for the effective 
average potential which we denote by dU kF /dt. Using the general formulae of [22] it can 
be computed without additional effort for the extended ansatz (3.3) where we keep the 
momentum dependence of Z^ tk and part of the momentum dependence of the Yukawa 
coupling with hk(q) = h(—q, q). With P F given in appendix D and setting for a moment 
Zip,k(q) — 1 one obtains: 



d_ 

di' 



^ = -2<- 1 ^/^E^(^)+^))' 1 • (3-23) 



Here m 2 a (q) are the N real nonnegative eigenvalues of the N x N matrix hk(q)h* k (q)ip^ip 
with momentum dependent Yukawa couplings defined by (A. 8). Here we have taken into 
account the N c colors of the quarks. For a given value of q we may use the identity (with 
d/dt acting only on P F , m 2 = diag(m 2 ) and tr taken in flavor space) 



d 



t^iPF + miy 1 = ^trln(P F W) = |lndet(P, + 



dt 



m 



— Indet [P F + h k (q) 



dP F 
~dt 



tr [P F + h k (q) 



This gives an expression for general 9? 



dt UkF 



d 

-22 



d d q dP F 
(2ir) d dt 



tr (P F + h k (q) 



(3.24) 



(3.25) 
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Using the particular configuration with N 
one finally obtains 



1 equal eigenvalues and the relation (3.13) 



d_ 

di 



U, 



kF 



d d q dP F 



J (27l) d 



(27T) d dt 



(N 



Pf + 



N 



hi 



+ 



hi 



(p-(JV-l)V^ 



(3.26) 



We emphasize that the Yukawa couplings in (3.3) conserve the axial ^7a(1) symmetry. 
The fermionic contribution to dU kF /dt is therefore independent of £. The fermionic wave 
function renormalization Z^^(q) is easily restored if we replace the function P F (q) in (3.26) 

by Z 2 k (q)P F (q) and note that the partial derivative ^(Z^ k P F ) only acts on the pieces 
related to the infrared cutoff R k - Within the truncation (3.3) the fermionic contribution 
to the evolution equation (3.25) is then exact. Eq. (3.26) gives the exact result for N = 2 
whereas for N > 2 one has an additional dependence on invariants Tj, i > 3, defined in 
appendix B, which can be extracted from (3.25). The derivatives of §fUkF with respect 
to p and r 2 can be written in a suggestive form as 

d d q 



dt kF 



-2 2 — 

N 



x — < 



d_ 

di 



(2n) d hk 
N-l 



(3.27) 



j_ 

N 



hi 



(P 



Z l,k P F + 



d dU, 



kF 



dt dr-y 



d 2A7 N - 1 
-2~ 2 N, 



d d q 

N (2nY 



hi 



(3.28) 



x — < 



d_ 

di 



Z%^Pf + ]v 



hk (p + V^) 



Z^k^F + Jj 



hi 



with the formal definition (cf. appendix D 

1) 1 dR k d 2 
01 



+ 



P F d [Z^ k r F ] d 



(3.29) 



Z^ k dt 8P Z^ k l + r F dt dP F 

We may now combine the bosonic and fermionic contributions to the running of the 
renormalized couplings. Here we restrict the discussion again to momentum independent 
Z<p, Z^ and (real) h, i.e. we replace similarly to (3.7) Z^^q) — > Z^k), Z^^ k (q) — > Z^(k) 
and h k (q) — > h(k) = h (k). For arbitrary d it is convenient to introduce dimensionless 
couplings in analogy to (2.8) — (2.10): 

h 2 

A 1 o 

(3.30) 



^1,2 — 



z~ 2 k a 



K 



e = 



k 2 ~ d PR 

1 -2 2 

k m 



z v k 2 - d Po 

Z- x k- 2 m 2 
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and to define the anomalous dimensions for the scalar field, r]^, and the fermion field, rj^, 
by 

- d \ 7 9 1 7 



(3.31) 



We also use dimensionless integrals 



n 



= -vl'k 



1 i2n-d 



d d q ( 1 dR k (q)' 



f d d q I 
J (2^p { 



Z v dt 



(P + wk' 



-(n+l) 



-k 



2n-d 



dxx 2 



(3.32) 



«( p+ »**)"" 

where the part ~ T] v l^(w) arises from the t-derivative acting on within Rk (cf. (1.2)) 
and 



1 = 2 d+1 J r f - 



(3.33) 



The "fermionic integrals" l^ d (w; rj^) = l^ d {w)—ri^l^ d {w) are defined analogously, with 
P replaced by Pp. Combining (3.14) — (3.16) with (3.27), (3.28) we obtain the evolution 
equations for the symmetric regime: 



de 
dt 

d\i 
~dt 

dM 
dt 



-(2 - Vv )e - 2v d {[2(N 2 + 1)X 1 + (N 2 - l)A 2 ]/f (e; V(p ) 

2iN c hH[ F)d M} (3.34) 

(d - 4 + 2 Vip )\i + 2v d { [2(N 2 + 4)Ai + (iV 2 - l)A 2 (2Ai + A 2 )]Z 2 d (e; % ) 



2 |* 



f^ 4 4 F)d (%)} 

(d - 4 + 2?^)A 2 + 2v d { 

2 i +1 §^f )d (w)} • 



(3.35) 



12AiA 2 + 2(N 2 - 3)A 2 Zg(e; 77^) 



(3.36) 



Similarly, the evolution equations for the SSB regime read 

_ = (2-rf-^) K + 2^{iV 2 /f(^) + 3/f(2A 1 / t ;^) 



+ (iV 2 -l) 



1 + 



A, 
Ai 



^(A 2 «;^)-2l7V c ^P(l^ ;?7V)) 



(3.37) 
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= (d-4 + 2 V ^\ 1 + 2v d {N 2 \jl d 2 (r ] ^ + 9\jl d 2 (2\ 1 K-,r ] ^ 

+ (N 2 - 1) [Ax + A 2 f l d 2 (X 2 K; r) v ) - 2^hH { 2 F)d (±h 2 K ] % )} (3.38) 

^ = {d-A + 2^)\ 2 + 2v d ^\ 2 2 l d 2 {^) + ^-{N 2 - 

iiV 2 A^f 5l (0, X 2 k; ^) + 3[A 2 + 4A 1 ]A 2 Z£ 1 (2A 1 k, A 2 «; Vv> ) (3.39) 

where we have defined 



l d M = l d (0; % ) , ^ d ( % ) = l^ d (0; th) . (3.40) 

No explicit dependence on k appears in this scaling form of the flow equations. In the 
limit 6, k — ► one recovers for both regimes to leading order in the coupling constants the 
known [16, 17] perturbative one-loop beta functions for Ai and A 2 . 

The system of flow equations (3.34)-(3.39) is the central piece of this work. For 
the quark-meson model (d = 4) we fix the initial conditions at k = k v (t = 0) with 
Ai(A^) = A 2 (/c ¥ ,) = 0, e(k ip ) = Z~ 1 (k (p )k~ 2 m . The solution of the evolution equations 
should then reveal the phenomenon of spontaneous chiral symmetry breaking for k — > 
(t — > — oo). More precisely, we expect for some scale k s > that the mass term vanishes, 
i.e. e(k s ) = 0. Subsequently, for k < k s we follow the evolution of the minimum of the 
potential using the system (3.37)-(3.39), with initial condition K,(k s ) = 0. (The couplings 
Ai and A 2 are continuous at k s .) For sufficiently small values of k the minimum of the 
potential will not move anymore, i.e. 

\\mZ v {k) = Z^ 

\im p R (k) = lim ZJk)p (k) = p R = Z^po 

k— >0 k^O 

N (3-41) 
= Na 2 R = NZy = -f 2 



\imK,(k) — > p R k 



-2 



Supplementing the flow equations for k, Ai and A 2 by the one for the Yukawa coupling 
h 2 and inserting the anomalous dimensions — these quantities will be computed in the 
next two sections — we can now study how the shape of the average potential changes 
as k is lowered. We have integrated the flow equations numerically for d — 4 and N c = 3 
from t i — corresponding to k = k v to tf = \n(m n / k^) . This endpoint of the numerical 
integration simulates the pion mass threshold which is neglected in our approximation 
of vanishing quark masses. We ignore here all dependence of the threshold functions on 
the anomalous dimensions rj v and rj^,. This approximations will be justified in section 7. 
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Figure 2: Flow of Ai (solid lines) and A 2 (dashed lines) with fc, for the Ul(2) x Ur(2) model with two 
sets of initial conditions at k v — 630 Me V with h 2 (k v ) = 300 and e = 0.01. 

We use first k v = 630 MeV, e = 0.01 and Z^(k ip ) = 1, /i 2 (£^) = Z~\k v ) = 300. For 
the initial values of Ai and A2 we employ two different sets of boundary conditions. One 
corresponds to the approximations used in [14] 

Ai(*v) = X 2 (k v ) = . (3.42) 

The other set is obtained in the large- N c limit of the Ul(N) x Ur(N) model [17] and 
reads 

Ai(M = |^ 2 (M 

4 (3-43) 

UK) = ^h\K). 

For the numerical investigations we will concentrate in the present work on iV = 2. 
In fig. 1 we have plotted the renormalized mass m as a function of k. Starting from 
a very large value m(k ip ) = 1091 MeV the mass rapidly decreases and reaches zero for 
k s ~ 450 MeV. For k < k s the minimum of the potential occurs for po > and we have to 
use the flow equations for the SSB regime. In fig. 1 we also show the /c-dependence of the 
location of the minimum, <r R {k) = (pR(k)fN)' 2 , and see how it stabilizes for k ^ 250 MeV. 
The final result for or is rather insensitive to the exact choice of the endpoint kf — m n . 
In fig. 2 we display the /c-dependence of the quartic couplings Ai and A2 for the two sets 
of boundary conditions. We observe that the result for k = kf does not depend strongly 
on the initial values. This is a first manifestation of the infrared stability mentioned in 
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the introduction. We will explain the origin of this behavior in more detail later (section 
7), since for an understanding we first need to discuss the anomalous dimensions and the 
running of h 2 . 



4 Scalar anomalous dimension 



A computation of f n in terms of the four-quark interaction at the scale requires the 
ratio Z^{$) / Z v {k v ). An evaluation of this ratio is the subject of this section. We will 
begin with the determination of the flow equation for the momentum dependent scalar 
wave function renormalization Z^ k {q) or, equivalently, the scalar anomalous dimension 
i] Vtk (q) = —Z~ 1 k {q)d t Z ip ^{q). For this purpose we have to consider a spatially varying 
scalar field configuration. We choose a nondiagonal distortion of the constant vacuum 
configuration (2.3): 



<p ab (x) = ip5 ab + 5<pe lQx + 5<p*e tQx S a6 = ip* ab (x) 



with 



S a 6 — 8 a l$b2 _ ^a2<5fcl 

or, in momentum space (with 5(qi,q2) = (27r) d <5(gi — 52)), 

^Pab(q) = (p5(q,0)5 ab + [5(p5(q,Q) + 5(p*5(q,-Q)]Y, ab 
= ip6(q,0)6 ab + A(g,Q)S a6 



(4.1) 
(4.2) 

(4.3) 



Expanding around ip at the potential minimum, we observe that Sip corresponds to an 
excitation of the massless charged pion tt ± . We keep the discussion here for general if. If 
we supplement the scalar configuration by a fermionic background 







(4.4) 



(2) 

T k is easily seen to be block-diagonal. It decays into matrices acting in scalar and 



fermion subspaces, T^l and T^, respectively. Hence, we can read off from (1.1) 



^(2) 



- Tr 



(Tpl + R Fk j 



-1 dR 



Fk 



dt 



(Q - 0)1 



-1 dR k 



dt 



(4.5) 



The right hand side may be evaluated by expanding the traces in powers of 5<f and 8ip* 
up to order 8ip5ip* and subtracting all Q-independent terms. The flow equations for the 
renormalization constant Z^{k) or equivalently r) v (k) are now determined as 

foZy(k) = \ua Q —Z V)k (Q) , r] v (k) = Hm^^Q) = -— hi Z v (k) . (4.6) 
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In order to evaluate the right hand side we split 

r (2) _ r (2) Ar (2) 

such that all Sip, Sip* dependence is contained in Ar^ 2 -*. We may then expand 



(4.7) 



Tr 



(rl 2) + i4) 



-i dRv 



dt 



Tr 



-i dRi 



+ Tr 
1 



iKrS + ^Arf)} 



Tr 



+ C(A 3 ) 



^ (4-8) 

and compute the right hand side of (4.5) from the first terms. Details of the calculation 
can be found in appendix E. Taking the limit Q 2 — > and neglecting all momentum 
dependence of the Yukawa coupling and wave function renormalizations we find for the 
SSB regime 



iV 2 



+ 22 



^A^mf^-U^) 



(4.9) 



~ AT 



and for the symmetric regime 



^ = 2f +2 ^iV c / i 2 mf )d (0;^). 



(4.10) 



Here we have defined the threshold functions 



{w 1 ,w 2 ;r] ip ) = m* un2 (w 1 ,w 2 )-ri v m* una (w 1 ,W2) 

= _k2(n 1+ n 2 -l)-d r dxx ld_i 

2 io <9t \ 



P(x) 



[P(x) + fc 2 wj ni [P(x) + k 2 w 2 ] n2 J 

(4.11) 



and 



m{ F)d (w; 7fy) = mf )d (w) - r]^m{ F)d (w) 



cfea; 2+1 — 



d ( d 1 + r F (x) 



where 



x = g 2 , P(x) = P(g) , P(x) 



<9t ycte Pf(x) + /c 2 w, 



— — Pm , 7— P = — — P 



(4.12) 



(4.13) 



From the anomalous dimension r] v (k) the wave function renormalization can be ob- 
tained by numerical integration of (3.31). We have plotted the result in fig. 3 for 
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Figure 3: Evolution of Z v with k, for the U L (2) x U R (2) model with k v = 630 MeV, e = 0.01, 
^i(k v ) = \2(k v ) = and two values h 2 (k v ) = 300 (solid line) and /i 2 (/c ¥ ,) = 10 4 (dashed line). 



k v = 630 MeV, e = 0.01, Ai(fc ¥ 



and two values of /i 2 (/c v ) = 300 and 



10 . We note that Z^ increases strongly in the symmetric regime for k > 450 MeV and 
stabilizes for low values of k. Again, the final value at k = kf does not depend very much 
on the initial conditions for h 2 . 



5 Evolution equation for the Yukawa coupling and 
fermion anomalous dimension 

To determine the evolution equation for the Yukawa coupling and the fermionic wave 
function renormalization constant we will turn to a field configuration (ip — <f*) 

(fiab(x) = (f5 ab 

C(x) = Ce"^ (5.1) 

Furthermore, we will approximate the momentum dependence of the Yukawa coupling by 
h)-(—q,q') ~ /ife( 2 -^ 2 -)- This amounts to neglecting its dependence on the external scalar 
momentum in the Yukawa vertex in (3.3). Accordingly, the renormalized Yukawa 
coupling is defined via 

h k (q) = 0- 2 Z^ k (q)Z;-i(O)h k (q) . (5.2) 
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The matrix of second functional derivatives of simplifies considerably for the above 
configuration. Omitting spinor indices one finds 

AT a/^ ; t = (Z+Mi +h k (q)Vl) (2n) d 5(q - q')5 b a . (5.3) 
dip {q)d^ b {q') v 7 

The derivation of the flow equations for h and follows similar lines as for the scalar 
anomalous dimension discussed in section 4. For details of the calculation we refer to 
appendix F. Neglecting the effects of the chiral anomaly (V = 0) as well as the momentum 
dependence of the wave function renormalizations and the Yukawa coupling we find in 
the limit Q — > 

JLtf = [d - A + 2 ^ + Vv }h 2 -^v d h^N% B)d (^h 2 ,e ]V ^ Vv ) 

- (N 2 - l)l[ F 1 B)d (^Kh 2 ,e + k\ 2 -,v^,V v ) -l[ F i B)d (-^Kh 2 ,e + 2K\ 1 -,r]Tp,r] ip )} 

(5.4) 

where 

CSV. ty) = CSVi, «*) - VKSV. «*) - VKVi> «*) 

= _I^.2(ni+n 2 )-d f°° dxX^' 1 — I - 1 

2 Jo dt\ [P F (x) + k 2 Wl ] n ^ [P(x) + k 2 w 2 } n2 J ' 

(5.5) 

Similarly, the fermionic anomalous dimension reads 

4 u d; 



* = Nd 



h 2 iN 2 m{ F 2 B)d {^-h 2 K, e; 77^, 77^) + m!™(^/i 2 /t, e + 2«Ai; 77^,, 77^) 

Ffld 1 (5 ' 6) 

+ (iV 2 - l)m{f )d (— h 2 K,e + K^ri^^^j 



with 



(FB)df 



l,2fn, +n ,-lW r , i^" f l+ffW P (^) 1 



_j^2(m+7i 2 -l)-<i / 

2 7o 9* I [Pf(z) + A; 2 wi]"i [P(x) + A; 2 w 2 ]" 2 J 

(5.7) 

In summary, the equations (4.9) and (5.6) constitute a linear system for 77^ and 77^ with 
solution 

A t (l + A 6 ) - A 3 A 4 

V<p = 



(1 + A 2 )(l + A 6 ) - A,A 5 

A^l + A^-A.A, 1 ' ' 

^ " (1 + A 2 )(1 + A 6 )-A 3 A 5 
where A ± , . . . , A 6 are defined by writing (4.9) and (5.6) in an obvious notation as 

= M- A 2 7} v - A 3 r}^ 

(5.9) 

77^ = A 4 - A^ - A 6 f]^ . 
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Figure 4: Dependence of h 2 on k, for the Ul(2) x Ur(2) model with two different initial values, 
/i 2 (fc v ) = 300 (solid line) and h 2 (k v ) = 10 4 (dashed line). We use k v = 630 MeV, e = 0.01 and 
Xilkp) = X 2 (k v ) = 0. 



We note that in four dimensions the integrals 



4f )4 (0,0)=mr^(0) 



m 



(FB)4 
1,2 



(0,0) = 1 



(5.10) 



are independent of the particular choice of the infrared cutoff. We therefore find in the 
limit of small masses k, e in both regimes to leading order in the coupling constants the 
known [16, 17] perturbative one-loop results for both anomalous dimensions: 



8tt 2 
N 



h 2 



(5.11) 



16tt 2 



h 2 



This in turn yields the correct perturbative one-loop result 



d_ 

dt 



h 2 = + 



N + N, 



8tt 2 



h 1 



(5.12) 



We observe that for large h 2 the running of the Yukawa coupling is very fast due to the 
term ~ h 4 in the flow equation (5.4). This explains why different large values of /i 2 (/c v ,) 
lead to very similar results for h 2 R = h 2 (0). We demonstrate this in fig. 4 where two 
different initial values /i 2 (/c ¥ ,) = 300 and /i 2 (/c ¥ ,) = 10 4 are compared. Fig. 5 shows the 
corresponding evolution of Z^(k). 
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Figure 5: Running of Z$ with k, for the Ul(2) x Ur{2) model with two different initial values h 2 (k v ) = 
300 (solid line), h 2 (k v ) = 10 4 (dashed line), and k v = 630 MeV, i a = 0.01, Xi(k v ) = \ 2 (k v ) = 0. 

6 The chiral anomaly 

So far we have considered the somewhat unrealistic limit V — > where the effects of the 
chiral anomaly are neglected. In view of the large value of vr (2.22) as compared to 
k v ~ 630 MeV, however, it appears that the opposite limit, V — > oo, should be closer to 
reality. For N = 2 it is straightforward to take the effects of the chiral anomaly in this 
limit into account. To see this we notice that the complex (2, 2) representation (p of the 
global symmetry group SU L (2) x SU R (2) ~ 0(4) decomposes into two irreducible real 
vector representations of 0(4) (cf. the discussion at the end of appendix B): 

V =\(a-tr ] ') + ^(a k + t7r k )r k . (6.1) 

By taking V — > oo while keeping m 2 — \u (or ~p 2 + ^T>) fixed, m a and m v ' diverge and the 
four corresponding mesons decouple. Hence, we are left with the real vector representation 
= (cr, 7Ti, 7r 2 , 7T 3 ) of 0(4). Its potential reads in the symmetric regime 

U k = \{m 2 - + l -M (0 a a ) 2 (6.2) 

and similarly in the SSB regime. We therefore end up with the 0(4) symmetric linear 
sigma model coupled to fermions. The flow equations in the symmetric regime are given 
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-(2 - Vip )e - 2v d {6Ai/f (e; Vv ) - 2* N c h 2 l[ F)d (^)} 
(d - 4 + 2r ]<p )X 1 + 2v d {l2A 2 Z 2 d (e; ^) - 2^ l N c h A lf ] \^)} 
(d - 4 + 2 W + ^)/i 2 - 4^/i 4 /i™(0, e; W , ^) ( 6 - 3 ) 



by [10, 23] 

de 

dt 
d\i 
~dt 
dtf 

dt 

fhp = 2l +2V -jN c h 2 mi F)d (0; VlP ) 

where e is defined here by e = Z~ 1 /c~ 2 (m 2 — |z7). For the SSB regime we find 

^ = (2-d-^) K + 2^j3/f(^) + 3^ 
9Ai 



dt 



dt 



= (d-A + 2r ]ip )X 1 

+ 2v d ^A 2 /^) + 9\*l*(2\ lK ; Vip ) - 2^ jV c fc 4 Zf )d (^ 2 «; ^) 

= {d-A + 2^ + r lv )h 2 (6.4) 

- 2^/i 4 |3Z^ B)d (^ 2 «, 0; - C^^^' 2A i K 5 ^)} 



^ = A^^AKXlml 2 (0,2X 1 K ] ^)+2lN c h 2 m{ F)d (h 2 K ]V ^ 



Vi> = 2— h |3mi >2 (-/i re, 0; 7fy, 77^) + mi j2 (-h re, 2Aire; 77^, 77^) 

The difference between the results of the 0(4) model (z7 — > 00) and the C/l(2) x C/r(2) 
model (17 = 0) can be taken as a measure for the uncertainty due to the rough treatment 
of the anomaly in the present work. Here the Ul{2) x C/r(2) model exhibits the effects 
of additional scalar degrees of freedom beyond the pions (and the cr-mode). Since the 
additional modes are relatively heavy, the 0(4) model should be closer to reality. The best 
model with N = 3, where the strange quark mass and the chiral anomaly are properly 
taken into account, is expected to deviate from the 0(4) model in the same direction as 
the Ul(2) x Ur(2) model. Since there are three light quark flavors in nature with masses 
smaller than k v one might naively expect the case N = 3 to correspond to the most 
realistic description of the real world. However, we have neglected quark masses and in 
particular the strange quark mass in this work. In the SUl(3) x SUr(3) x C/y(l) model 
the four if-mesons will therefore appear as massless Goldstone degrees of freedom which 
will unnaturally drive the evolution of all parameters even at scales much lower than 
their physical masses of approximately 500 MeV. The same holds, of course, for the three 
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pions. The effects of their physical masses m w ~ 140 MeV can, however, be mimicked by 
stopping the running for kf — m n . We therefore expect the case N = 2 to yield more 
realistic results than N = 3 as long as the strange quark mass is neglected. In addition, 
we note that for N = 3 the scalar self coupling Ai turns negative for positive but small 
values of m 2 . This happens despite the fact that Ai has acquired first a large positive 
value due to the strong initial Yukawa coupling. The cause is a large value of A2 which 
can drive Ai negative when the scalar loop contributions to the running of Ai become 
numerically important around the scale k s . We interpret this quartic instability of the 
truncation (2.7) of the effective potential as a signal for a first order phase transition in 
the mass parameter even within the Ul(3) x C/r(3) model without chiral anomaly. A 
proper treatment of the case N = 3 therefore requires a more general truncation of the 
effective potential [24]. We will leave this problem for future work. However, we would 
like to point out that the inclusion of the cubic (for N = 3) £-term into the potential 
will change the phase transition to first order anyway. One expects that the vacuum 
expectation value of tp jumps discontinuously to a finite value already for a large value of 
e such that the scalar contributions to the evolution of Ai never become strong enough to 
turn it negative. 

The difference between the 0(4) model and the U L (2) x U R (2) model is exemplified 
in fig. 6. There we show the "prediction" of f w as a function of the initial scale k v . The 
difference between the two models is not very large. We observe for not too small values 
of ktp a linear behavior f n ~ k v . For fixed initial values of the dimensionless parameters 
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like e = m 2 Z^/k 2 h 2 etc., this proportionality follows on pure dimensional grounds if no 
other mass scale is present. The bending of the curves for small k v is therefore purely a 
consequence of the additional infrared cutoff kf — m w which obviously plays a negligible 
role for a realistic size of k^. Looking at fig. 6 the reader may prematurely conclude 
that the predictive power for /„. is severely limited by the arbitrariness of the choice of 
the scale k v from where on the meson description is used. We should emphasize that 
for a full treatment along the lines of [14, 15] this is actually not the case: If one lowers 
the transition scale k v more fluctuations are included in the momentum range q 2 > k 2 
where a quark-gluon picture is used. As a result, the pole like structure in the effective 
four-quark interaction becomes stronger and m 2 therefore decreases. This is the same 
effect as found in the quark-meson picture used for the fluctuations with q 2 < k 2 : the 
mass term m 2 decreases with smaller fcasa result of the Yukawa coupling to the quarks. 
In the limit where the pole like structure dominates the evolution of G(0) = in 
the quark gluon picture the running of fn 2 is identical in both pictures for k larger or 
smaller than k v . The initial value fn 2 (k v ) as a function of k 9 follows therefore the same 
renormalization group trajectory as given by the flow equation (3.14). In this ideal case 
the choice of the transition scale k v does not affect the results, since the initial conditions 
move on trajectories of constant physics. In practice, this ideal scenario will often not 
be fully realized, since different types of fluctuations are included in the quark-gluon and 
the quark-meson description. The dependence of the results on k v can then be used as a 
quantitative check of the reliability of the employed truncations for the effective action. 



7 Infrared stability and predictive power for f n 

Comparing the results for /„. from fig. 6 with the experimental value /„■ = 93 MeV we find 
a surprisingly good agreement for k v = 630 MeV as infered from ref. [14]. The question 
arises to what extent this result depends on the particular choice of initial values at the 
scale kp. In principle, the values of the parameters of the quark-meson system at the 
scale ktp can be computed from QCD [14, 15]. In practice, however, many quantities will 
not be available with high accuracy, since one has to deal with a problem involving strong 
interactions. If f n would depend very sensitively on such quantities, a computation of f n 
with satisfactory precision would be extremely difficult. We will argue in this section that 
for small enough the opposite situation is realized. In this event the prediction for f w 
turns out to be almost independent of the initial values of many couplings. The reason is 
that a small Z v corresponds to a strong Yukawa interaction. The large value of h induces 
then a very fast running of almost all couplings towards values which are determined by an 
infrared attractive behavior. More precisely, the ratios \\/h 2 and A2//1 2 are determined by 
infrared fixed points of the type first found in the electroweak standard model [25]. This 
explains the insensitivity with respect to the initial values as demonstrated in fig. 2. The 
Yukawa coupling itself is also strongly renormalized and predicted to be in the vicinity 
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of the upper bound of the relevant infrared interval 5 [26]. Here the upper bound of the 
infrared interval is essentially determined by the scale k s where spontaneous symmetry 
breaking sets in and the effective quark masses constitute an infrared cutoff. This can 
clearly be observed in fig. 4 where also the insensitivity with respect to the initial value 
h{ktp) becomes apparent. The only relevant parameter will turn out to be the ratio 



e(fcy) _ m 2 ( y k !p )Zl(k v ) 
l> 2 (k.) ~ 



U2 



(7.1) 



This value determines k s and f n as well as all other couplings at the scale kf. 

For small enough Z v one starts in a regime where e = m 2 /k 2 = m 2 Z~ l k~ 2 is large. 
For large e the scalar fluctuations are suppressed by inverse powers of e appearing in 
the threshold functions. Then the scalar fluctuations can be neglected and only quark 
fluctuations drive the flow of the couplings. This is the approximation used in ref. [14] 
which remains valid as long ase > 1. The Yukawa coupling h is normalized according to 
(3.7) as h(0) = /ifc (0, 0) = 1. Consequently, the initial value for the renormalized Yukawa 
coupling is 

We use here a normalization of the fermion kinetic term such that Z^(k^) = 1. For small 
Z v {k v ) we therefore start with a strong Yukawa coupling. In the limit e ^> 1 the flow 
equations simplify considerably. If we define e = e/h 2 and A« = \i/h 2 we find 



de 
dt 

~dt 

8X2 
dt 

dt 



-(d-2)i+2l +1 v d N c 



= 2i +l v d N c K 2 



d 



Ai 



1 

N 



2i +2 v d N c V 



i X2 ~ k 



(7.3) 



(d-4)h 2 + 2-* +2 ^N c h 4 



In the following we will specialize to the case d — 4. 6 

As a first observation we notice that the Ai-A 2 system exhibits an infrared fixed point 
given by 

Ai* = \^2* = ^ ■ (7.4) 

This fixed point corresponds exactly to the large- N c estimate of [17]. The explicit solution 

5 The infrared fixed point for h is /i* = if no infrared cutoff is present. Due to a finite amount of 
running from k v to k s this translates into an infrared interval. 
6 The system (7.3) remains solvable for general d. 
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of the differential equations (7.3) reads 



m 

~X 1 (t) 

h\t) 



Av A N c 



1-e 



-it 



+ e e 



-2t 



A 10 



- 8§v 4 h 2 t 
1 - AN c v 4 h 2 t 



A20 



- lQ§v A Bt 



N 



(7.5) 



AN c v 4 hlt 
h 2 



I - AN c v 4 h 2 t 



with eo = e(t = 0), etc. The system crosses into the SSB regime when e(t s ) 
corresponding to a scale 



ts 



I In 

2 



4iV> 4 



k 2 



1 



e fc, 



= 0, 

(7.6) 



Here we have used v 4 = (327T 2 )" 1 and N c = 3 in the last expression. Around the scale 
k s our approximation (large e) breaks down. Nevertheless, it becomes apparent already 
at this stage that values of eo substantially larger than 0.04 are incompatible with chiral 
symmetry breaking, since e would remain positive for all k in this case. For large eo 
the effect of the quark fluctuations is simply not strong enough to turn the scalar mass 
term negative. We notice that e < 1 for m(k v ) ~ O(fOOMeV) and Z^k^) ^ f. If we 
furthermore assume — 4iY" c f 4 /iQt s ^> I or 



e >4iV>4 I-exp{- 



2N c v 4 h 2S 



(7.7) 



we find for not too large \i /hl, A 2 o/^o 



Ai(t s ) ^ Ai* , X 2 (t s ) ~ \ 2 * , h 2 (t s ) ~ 



2N c v A In 



4iV c v 4 



n -1 



4iV> 4 - g 



(7.8) 



This result may be interpreted as follows: Even though (7.6) should only give an approx- 
imate estimate for the entering point into the SSB regime this is sufficient to imply that 
\i and A 2 approximately reach their fixed point long before e goes through zero provided 
(7.7) is fulfilled. Furthermore h 2 becomes asymptotically independent of and is only 
a function of eo. AVe therefore conclude that for small Z lfi (k ip ) the system is governed by 
an infrared fixed point of Ai and A 2 in the symmetric regime and looses almost all its 
information on the initial values at k = k v . Infrared quantities of the SSB regime like 
or meson masses will therefore merely depend on eo for a given k v . It is straightforward 
to see that the same analysis holds for the 0(4) model discussed in section 6. The ap- 
proximate flow equations for e, Ai and h 2 are the same as (7.3). The difference between 
the Ul(2) x Ur{2) and the 0(4) model arises principally from the behavior of the running 
couplings around the scale k s . 
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Figure 7: The 0(4) model pion decay constant / T as a function of h = h(k ip ) for k v = 630 MeV, 
eo = 0.02 and Xi(k v ) = (solid line) as well as Xi(k ip ) — (dashed line). 



For larger values of Z ip (k ip ) or smaller values of h{k v ) the attraction of the infrared 
fixed points becomes weaker. As a consequence, the dependence of f n on the initial values 
of the couplings becomes more important as demonstrated in fig. 7. We conclude that 
for h{ktp) substantially smaller than ten it will become more and more difficult to obtain 
an accurate prediction for /„.. On the other hand, fig. 7 clearly shows the approximate 
independence of f n on h 2 (k v ) or \i{k v ) if h 2 (k,p) exceeds 300. 

An additional aspect of strong Yukawa couplings concerns the error in /„- due to the 
truncations of the quark-meson effective action. The effects of truncations in the scalar 
sector are diminished by the fact that scalar fluctuations are subdominant in the region of 
very large Yukawa couplings. A similar argument justifies the approximation of neglecting 
the terms proportional to the anomalous dimensions in the threshold functions. For the 
m-type functions this approximation is valid, since ^2,^3,^5,^6 «C 1 in (5.8), (5.9). 
For /-type functions one might be worried that i] v and 77^, are large for the initial part of 
the running in the symmetric regime due to large values of h 2 . However, because of the 
large values of e in this range of scales the contributions from the /-type functions can be 
neglected altogether. For smaller values of e, i.e. for scales closer to k s , the anomalous 
dimensions are expected to be already small (as indicated by figs. (3) and (4)). 
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Figure 8: The pion decay constant /„■ as a function of e for k v = 630 MeV, Xi(k lp ) = \2(k v ) — and 
h 2 {k v ) = 300 (solid line) as well as h 2 (k v ) = 10 4 (dashed line). 

8 Computation of f n for strong Yukawa coupling 

In the last section we have shown that for large enough Yukawa couplings, say h 2 ^^) > 
200, the value of f n only depends on the parameter e(k v ) = e for given k v . We demon- 
strate this quantitatively in fig. 8 where we plot f n as a function of eo, for both the 0(4) 
and the Ul(2) x Ur(2) model as well as for two different large initial values of h 2 . The 
"prediction" for f n is rather insensitive to /i 2 (/^) for moderate values of eo. 7 For the 0(4) 
model we also observe an extended plateau where f w is not very sensitive to e either. For 
this plateau the value of f n comes out between 80 and 100 MeV which fits very well with 
the experimental value of 93 MeV. Also the renormalized Yukawa coupling Hr = h(kf), 
or, equivalently, the constituent quark mass m q = \h^f^ depends essentially only on eo- 
We show this in fig. 9, again for two different large values of /i 2 (/c (/p ). Since both Hr and /„. 
are functions of only one parameter e , there arises a quantitative relation between those 
two quantities. Consider first the 0(4) model: For a constituent quark mass of 300 MeV 
or h,R ~ 6.5 we read off from fig. 9 that e ~ 0.02. Inserting this into the plot of fig. 8 
one obtains 

f n ~92MeV. (8.1) 

7 Thc decrease of f„ for h 2 {k v ) = 10 4 observed for unnaturally small values of eo results simply from 
the fact that the system enters almost immediately into the SSB regime, having no "time" for Z v to grow 
much beyond Z lp (k v ). 
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Figure 9: The renormalized Yukawa coupling as a function of eo for k v — 630 MeV, Xi(k v ) = 
\ 2 (k v ) = and h 2 {k v ) = 300 (solid lines), h 2 (k v ) = 10 4 (dashed lines). 

This value was obtained for /i 2 (A^) = 300 and \i(k lf ) = \2(k ip ) = 0, but it turns out 
to be not very different for h 2 ^^) = f0 4 or different initial values of Ai and A 2 . A 
similar procedure gives for the U L (2) x U R (2) model a value e ~ 0.008 and in turn 
fx ~ 126 MeV. Repeating this procedure for k v = 700 MeV we can infer from table 1 a 
guess of the uncertainty in /„- that can be expected within the quark-meson model for 
large Yukawa couplings, as represented in the table by the values h 2 (k,p) = 300 and 10 4 . 
On the other hand, a given value of m q also implies a minimal value h m i n (k v ) such that 
the evolution of h(k) can reach the value h R = 2m q / f w at all. We assume here that the 
result for m 2 (k v ) of ref. [14] should not be off by more than a factor of four. We can 
therefore conclude that for k v in the range (630 — 700) MeV there exists a lower bound 
on eo, i.e. e"n ^ 0.01. This in turn amounts for the 0(4) model to /i m i n (^) — 6.2, 8.0, 10.3 
for m q = 250, 300, 350 MeV, respectively. The corresponding values in table 1 give an 
estimate for the maximal deviation of /„- from its value for strong Yukawa coupling. For 
the Ul(2) x Ur(2) model a value h 2 (k,p) = 300 is already near the lower limit of what is 
compatible with realistic values for m q and e (cf. figure 9). 

We can also invert these relations and look for the optimum value of eo and /i 2 (/c </3 ) for 
fixed U = 93 MeV and m q = 300 MeV. One obtains for the 0(4) model and k v = 630 MeV 

e ~ 0.02 

(8-2) 

h 2 {k v ) ~ 280. 

Within the simple QCD inspired model of ref. [14] the transition scale was found as 
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m q 
McV 


300 


350 


h' 2 (k in ) 


64 


300 


10 4 


106 


300 


10 4 
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McV 
















0(4) 


630 




143.9 


91.7 


83.5 


124.6 


99.9 


91.0 


700 




159.7 


101.5 


92.5 


138.2 


110.7 


100.7 


U L (2) x 17,2(2) 


630 






125.7 


138.3 






138.9 


700 






139.6 


153.5 






154.2 



Table 1: in MeV for various initial conditions at k v and two values of the constituent quark mass. 



hp = 630 MeV and the mass term at k v gave e"o = 0.036. Comparing with fig. 8 we find 
that this value of e is actually too large to induce spontaneous symmetry breaking if the 
meson fluctuations are taken into account. Given the simplified character of the model 
considered in [14], however, we find the agreement with the order of magnitude of (8.2) 
very encouraging. On the other hand, the estimate of Z v (k v ) ~ 0.85 appears to be very 
inaccurate for the model of [14] and far away from the small values of Z v (k v ) for which 
the values (8.1), (8.2) were obtained. 

We finally compute the chiral condensate from (2.13) with h(k,p) = 1 as 



e Z v 2 {m % )Z i> (k v )f n k 



(8.3) 



Extracting Z^im^) from section 4 (fig. 3) we can use this value for a check of the self- 
consistency of our scenario. We observe that three different small quantities, eo, 2,,(m,) 
and fn/ktp enter here, and it is by far not trivial that a reasonable value of the chiral 
condensate can be obtained. We normalize the condensate at k v with Z^ik,^) = 1. The 
result for (ipip) as a function of e is plotted in fig. 10. For k v = 630 MeV, e = 0.02 



and h 2 (k v ) = 300 we obtain in the 0(4) model 



~ 163 MeV . 



(8.4) 



We may compare this value with a typical value infered from chiral perturbation theory 
[3] 



CPT 



1 GeV) ~ (225 ± 25) MeV 



.5) 



This value can be scaled down to k u 



630 MeV by exploiting k-^ (^) CpT (k)m q (k) = 0. 



We use here the three-loop /3-function of QCD (MS scheme) for the running quark mass 
[27] 



k d 
m q {k) dk 



m q (k) = 



a s 

7i \-l2 

7T 



+ 73 



IT 



(8.6) 
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Figure 10: The quark condensate as a function of eo for k v = 630 MeV, \i(k v ) — A 2 (fc v ) = and 
h 2 (k v ) = 300 (solid line) as well as h 2 (k v ) = 10 4 (dashed line). 



with 



7i = \c F \ 72 = [97iV c + 9C F - ION] 
Z 4o 



73 32 



11413 2 129 
108 c 



(8.7) 
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(23 - 24C(3))iVC F -^iV 2 



and C F = %±. This yields 



for = 630 MeV and therefore 



m 9 (lGeV) 



1.72 



^) CPT 3 (M^ (188 ±21) MeV. 



.9) 



This is in satisfactory agreement with our estimate. 



9 Conclusions 



We have presented here an effective quark-meson model which is supposed to describe 
the strong interaction dynamics between (constituent) quarks, scalar and pseudo-scalar 
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mesons at momentum scales smaller than k v ~ 630 MeV. The effective average action 
T k for this model depends on a scale k which plays the role of an infrared cutoff. The 
scale dependence of the average action obeys an exact nonperturbative evolution equation. 
Using a truncation for the general form of this results in approximate flow equations 
for the meson potential and kinetic term as well as the quark kinetic term and the quark- 
meson Yukawa coupling h. The initial values of these parameters at the scale k v can, 
in principle, be computed [14, 15] from evolution equations for QCD which are valid 
for scales larger than k^. Following the flow equations from k^ to lower scales k one 
recovers for k = the effective action, i.e. the generating functional for the 1PI Green 
functions for the mesons. In particular, the standard nonlinear cr-model framework of 
chiral perturbation theory obtains if the expectation value of the meson field is kept fixed 
at a nonvanishing vacuum expectation value. 

We have solved the flow equations numerically and observe how the minimum of 
the meson potential turns from <t = at high scales to a nonzero value | cr 1 > for 
small k. The nonvanishing expectation value <r indicates spontaneous chiral symmetry 
breaking. Our numerical solution allows us to compute the mass scales characteristic 
for chiral symmetry breaking, i.e. the pion decay constant the constituent quark 
mass m q and the chiral condensate (i/ji/jJ^. These quantities are computed as functions 
of the initial values for the scalar mass term m 2 ^) and wave function renormalization 
Z^kp). The latter is related to the value of the renormalized Yukawa coupling at k v by 
h 2 (k<p) = Z~ 1 (k ip ) once we normalize the quark wave function with Z^k^) = 1. 

The large ratio between the constituent quark mass m q and f n necessitates a large 
value of the renormalized Yukawa coupling h 2 R = h 2 (k = 0) according to 

o 4m 2 , 

J 7T 

Since h 2 (k) decreases rapidly with decreasing k we conclude that at the scale k v we have 
to deal with a strong Yukawa coupling. Our investigation therefore concentrates on large 
initial values /i 2 (/c ip ) <^ 200. The most crucial observation of the present work is that 
strong Yukawa couplings imply a very fast running of almost all couplings towards values 
determined by infrared fixed points or corresponding infrared intervals. The quark-meson 
model for small scales k therefore looses its memory of the exact initial values of most of 
the couplings at the scale k^. In consequence, strong Yukawa couplings greatly enhance 
the chances for a reliable estimate of f w \ We find that f w as well as m q and (^pip^ depend 
only on one "relevant" initial value, namely eo ~ m 2 ^)/^ 2 . The value of h R may then 
be used to fix the required initial value of eo and therefore to determine /„- and (ipip 



o 

For the simplified 0(4) model discussed in the present paper we obtain /„. ~ 92 MeV 



and 



iptp) 3 — 163 MeV. The comparison with the experimental result = 93 MeV 



and the estimate from chiral perturbation theory yP^/^^ (k<p) = (188 ± 21) MeV is very 
encouraging! 

Despite the success of the simplified computation of the present work our results have 
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partly the character of a feasibility study. Several steps should considerably improve the 
accuracy of our computation of f n . First, one should consider the case N = 3 with 
a realistic value of V for the chiral anomaly. This investigation is already prepared in 
the appendices of the present paper where the central results are exhibited for arbitrary 
N and V. Second, the explicit chiral symmetry breaking due to current quark masses, 
in particular the strange quark mass, should be included. Third, the dependence of the 
Yukawa coupling and the quark kinetic term on the quark momentum can be incorporated. 
This will provide information on the quark wave function inside the mesons [14]. 

Furthermore, we have not included residual gluon effects for k < k v in the present 
work. This would not be necessary if the gluons could completely be integrated out for the 
determination of the effective quark-meson action at the scale k v . Since the latter seems 
to be a quite difficult task, one may rather use the proposal of ref. [15] and integrate out 
only the gluons with momenta q 2 > k 2 The remaining gluon fluctuations with momenta 
q 2 < k 2 give then additional contributions to the flow equations in the quark-meson 
model. In the context of a computation of /„- the most important effect seems to be a 
residual gluonic contribution to rj^ and the effective quark meson vertex and therefore to 
the running of the Yukawa coupling. This can be taken into account by generalizing the 
average action to include gauge bosons [28] and using the formalism of [15] to integrate 
them out consecutively. A great part of the gluon fluctuation effects is already included 
in the contributions from the effective four-fermion interactions for k > k v or from the 
quark-meson interactions considered in the present paper for k < k^. The corrections 
from residual gluon fluctuations can be found from the explicit formulae in [15]. Since in 
addition the confinement scale is below the constituent quark mass one may hope that 
the complicated effects of gluon condensation do not have a very important influence on 
the determination of 

Finally, a computation of /„- within QCD, i.e. as a function of a s (l GeV) and the quark 
masses, necessitates a reliable computation of e and /i 2 (/c v ,) within the QCD framework for 
the effective average action proposed in [15]. Thus, the road to an analytical computation 
of fit{a s {l GeV)) is still long. Our results should encourage to go it. 
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Appendices 



A Linear a— model with fermions 

In this appendix we describe the Yukawa couplings of quarks to mesons in the language 
of the linear <7-model in a Euclidean formulation. In Euclidean space a Weyl fermion is 
described by two independent Grassmann variables ip a and rp & . (We employ the notation 
of ref. [22].) We therefore describe the left-handed quarks and their antiparticles, the 
right-handed antiquarks, by spinor fields ipia and ip LA , the whereas right-handed quarks 
and left-handed antiquarks are contained in ip R and ip R , respectively. In this language a 
Dirac spinor is composed from left- and right-handed Weyl spinors: 



The tensors contained in the fermion bilinears are described in the usual way. For example 
a vector reads 

u m \ ( IpL \ 

? m )\1> R ) (A.2) 

The chiral SUl(N) x SU r (N) transformations (with N being the number of quark flavors) 
act independently on the left- and right-handed quarks: 



U R 



1> 




(A.3) 



u\ 



L J 

and the Ur(1) symmetry corresponding to conserved baryon number acts as 

i> -> expl+— jip 

( iM- (A ' 4) 

ip -> exp ( j ip . 

Here Ul and Ur are SU(N) matrices with U\Ul = UrIIr = 1. The vector (A.2) and 
therefore also the kinetic term for the quarks is invariant with respect to these transfor- 
mations. On the other hand, no mass term is allowed by the symmetry (A.3). 

Let us now introduce a complex scalar field ip in the (N, N) representation of SUl(N) x 
SUr(N) carrying no baryon number. The appropriate flavor transformations are 

* ~* U *M (A . 5) 
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where ip is represented by an N x N matrix. Its iV 2 complex components should describe 
mesons or quark-antiquark bilinears. The most general Yukawa coupling between mesons 
and quarks which is consistent with the SU L (N) x SU R (N) x U B (1) flavor symmetry and 
Osterwalder-Schrader positivity reads 8 

In general, the Yukawa couplings may contain form factors. In a momentum space repre- 
sentation one has 

S Yuk = nJ2 {h(-q, q')^ R (qMq ~ 4)M<J) - h*(~q, q')^ L {q)^{q' ~ q)M4)} ( A -7) 
g,<i' 

where we work here on a torus with volume Q such that momentum integrals are discrete 
sums. In particular, this yields the couplings to a constant scalar fields ((p(p) = <p6(p, 0)) 

S Yuk = Sl^(q) ( Kq ^ _ r ° {q)(ft ) m (A.8) 

with h(q) = h(—q,q). We adopt conventions with real Yukawa coupling h. By a suitable 
chiral transformation of fermionic phases this is equivalent to the convention with purely 
imaginary h which would be closer to the most commonly used phase convention in a 
Minkowski space formulation. 

By virtue of the SUl(N) x SUr(N) symmetry we can express the scalar potential 
Uk((fi) as a function of a small number of invariants. The number of independent invariants 
depends on A" and we want to classify them. For this purpose it is useful to consider in 
addition the axial Ua(X) transformation 

i) L -> exp(+i'd A )'ip L , i\) L ->■ ex.p(-i'& A )ip L 
ip R exp(-n? A )^fl , i> R -> exp(+i<d A )ip R (A. 9) 

ip — > exp (— 2id A ) 99 , — > exp ( y +2i r d A ) . 

(We observe that the Yukawa coupling (A.6) is invariant under U A (1) transformations.) 
With the help of SUl(N) x SU r (N) x U a (1) transformations we can bring (p into a 
standard diagonal form with real nonnegative eigenvalues m a . This implies that the 
SUl(N) x SU R (N) x U A (1) invariants can always be expressed in terms of the A" eigen- 
values rh 2 a . There are therefore exactly A" independent invariants. Without the U A (l) 
symmetry the standard form of ip is exp (i$) diag (jn a ) and the additional phase corre- 
sponds to an additional invariant with respect to SU L (N) x SU R (N) which is not invariant 

8 In the conventions of rcf. [22] we take e — 1 and use t)l = tpL, VR = V'.R; ^lr — h R L = h, 
h LR = h RL = —h*, hjii = liLR = h' LR = h' RL = 0. In this language the Yukawa couplings hij are 

2N x 2N matrices of the type ^ ] or | ^ | . 

V h J \ h 
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under Ua(X) transformations. If the symmetry <p — > </r is preserved we may choose for 
this invariant 

f = det up + det <p f . (A. 10) 

Consider now the A" SU L (N) x SU R (N) x C/a(1) invariants corresponding to the A" eigen- 
values m 2 a . The first of these invariants is simply the trace 

p = tr(^V) =E™«- (A.11) 

a 

For N > 2 this is the only invariant quadratic in the field ip — there is only one singlet 
in the direct product (N, N) <g> (N, N). Only for A^ = 2 the invariant £ is also quadratic 
in ip. It corresponds to the singlet in (N, N) (g) (N, N). The traceless matrix 

= «^V_l p (A.12) 

can be used to construct higher invariants, i.e. 

f; = tr0\ (A. 13) 

The invariants fj can easily be expressed in terms of the m 2 , using identities of the type 

tr 2 = tr (ip V) 2 - 1 (tr y. V) 2 = £ f% - p l (A. 14) 

tr0 3 = E^a " ^Ptr0 2 - ^p 3 . (A.15) 

For N = 2 the invariant 73 can be expressed as a function of p and fi and similarly for 
all higher %. More generally, we only need to consider the A" invariants p, ?2, . . . , fjv and 
find the most general SUl(N) x SUr(N) symmetric scalar potential 

U k (f) = U k (p,f 2 ,...,f N ,0 ■ (A.16) 



B Scalar mass spectrum 

In this appendix we compute the scalar mass spectrum in dependence on the background 
field for a general form of the potential. We start with an SUl(N) x SUr(N) x Ua(1) 
invariant potential. Without loss of generality we can then consider configurations of the 
type 

P>ab = VaSab = m a 8 ah . (B.l) 

We parameterize the potential C4(<p, <pt) in terms of the invariants p, Tf 

U k (<p,<p*) = U k (p,n) (B.2) 

p = tr(<pV) 

A" N / f 1 

N 2 1 
7-3 = (iV-l)(iV-2) f3 + r2 (R3) 
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whose relation to the eigenvalues m 2 a can be found in appendix A. (If needed, suitable 
definitions of t { in terms of f { have to be chosen for % > 4.) One obtains for the second 
derivatives of the potential 

62Uk U' k 5l5i + 2U'^ c 5 ab 6 cd + E ^ ^ 



Red 



U k 6Jb + E , ( B - 5 ) 



; Cft = 0. (B.6) 

Here </2 a ft = ^ (y?R a 6 + i<Piab) an d primes denote partial derivatives with respect to p. 
We observe that Sri/Sip^ vanishes for the configuration (B.l). In order to gain a better 
understanding of the mass matrix (B.4) — (B.6) we briefly discuss a few special cases: 

For the origin at tp a — we are in the symmetric regime and the mass matrix has 
2N 2 real eigenvalues U' k (0, . . . ,0). Here we use the fact that the invariants t { are at least 
quart ic in ip. 

If the potential is independent of Tj it exhibits an enhanced symmetry SO(2N 2 ) instead 
of SUl(N) x SUr(N) x C/a(1)- In case of spontaneous symmetry breaking the minimum 
of the potential occurs at U' k (po) = 0. The N 2 x iV 2 matrix ip a l fc^ab^ cd has exactly one 
eigenvalue p whereas all other eigenvalues vanish. Together with the massless fields (pi 
one therefore finds for p = po the expected 2N 2 — 1 massless Goldstone bosons. The radial 
excitation has mass squared 2U k (po)po- 

We next include the dependence on the invariant r 2 = j^zi^ = jjzi tr (vV) ~~ Jf^iP 2 - 
For the configuration (B.l) one has 

2V26 cd (-^- V l--^—pp c ) (B.7) 



Red 



i5:5t(—y a+ rf)-—p) 



2N 4 

+ J^Va<PbS d a S C b - J^—^^eSabS Cd (B.8) 

JWk7ed = ^ (^T ( ^ + ^ ' iv^T P ) - • (R9) 

The mass matrix for the fields (fj is given by (B.5) and easily evaluated if possible con- 
tributions from Tj, i > 2 are neglected. The N fields pi aa do not mix with (fi c d, c ^ d. 
The corresponding N eigenvalues of the mass matrix are 
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For a ^ b the fields (pi a b and ipiba mix, but decouple from ipjd if c 7^ a or 6 or if d 7^ a or 
6. There are 7V(iV — l)/2 eigenvalues 



2 g£4 
TV - 1 dT 2 

and iV(iV — l)/2 eigenvalues 

.2 2 dU h 



N (y 2 a + if 2 b + <p a <p b ) -p]+U' k 



(B.11) 



K + a6 )' 



iV — 1 <9r 2 



iV(^ + ^-^)-p]+^- (B.12 



The discussion for <pr is similar. For a ^ b only and (finba m i x an d one finds 

N(N — l)/2 eigenvalues 



( M Rab) 2 = J^ZY^- [ N + ^ + *W / ; 



(B.13) 



and N(N — l)/2 eigenvalues 

,2 2 <9f4 



N - 1 <9r 2 

The mass matrix for the TV fields <fR aa reads 



jV^ + y^-y^-pJ+f^. (B.14) 



M Rac 



u f k + 



2 8U k 
N - 1 dro 



+ 



4 <9 2 £4 
(TV - l) 2 (9r 2 ) 



<9r 2 



(B.15) 



(9r 2 



Its eigenvalues M|, ac have, in general, no particularly simple form. 

If we specialize to N equal values ip 2 a = jjp the matrix M\ simplifies considerably: 



M 2 Rac =lUl + 



4p dUk 
K N- 1 0t 2/ 

This matrix has (N — 1) eigenvalues 



5ae + 2U' k ' 



4 0tV 



N -ldr 2 
4p 017* 



and one eigenvalue 
For this special case one finds 



k ' N - 1 dr 2 
(M*y = U' k + 2U k > p. 

(mf*) 2 = ( M Rab) 2 = K) : 



(B.16) 

(B.17) 
(B.18) 
(B.19) 
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and 

M| a = (M+ & ) =(M i?ab ) = E/£. (B.20) 

For C/j(. = the iV 2 massless Goldstone bosons (B.20) correspond to the symmetry break- 
ing SU L (N) x SU R (N) x U A (1) -> S77(JV) where the unbroken S77(jV) is the diagonal 
subgroup of SUl(N) x SUr(N). In addition we have A^ 2 — 1 massive scalars (B.19) in the 
adjoint representation of SU(N). For U' k (po) = their mass terms are positive provided 
dUk/dT 2 > 0. Finally there is a singlet with mass term 2L^'(po)Po- 

Another interesting special case occurs for <p 2 = p8 a \ which corresponds to the symme- 
try breaking SU L (N) x SU R (N) x U A {1) -> SU L (N - 1) x SU R (N—1) x 17(1) x 17(1). The 
eigenvalues M| a contain A^ — 1 values — and one value + 2p^. The eigen- 

values (M+ b ) 2 , (Mj ab )\ (M+ ab ) 2 and (m^ 6 ) 2 decompose each into (N - 1)(JV - 2)/2 
values U' k - and JV - 1 values E/£ + 2p§^. Finally M| ac has N - 1 eigenvalues 

U'k-wh^ and one eigenvalue ^ + 2p^ + 6pg Wg + Sp 3 Jg^. For ^ + 2p^ = 
we observe the expected AN — 3 massless Goldstone bosons. 

Next we should take into account that the axial Ua(1) symmetry is broken due to 
anomalies. The most general potential can therefore also depend on 

f = det tp + det <p f . (B.21) 

We observe that £ violates {/a(1) but is invariant under SUl{N) x SUr(N). The second 
possible C/a(1) violating SUl(N) x SUr(N) invariant, = % (det <p — det <p^, violates the 
discrete symmetry ip — > </r and hence the CP invariance of the model and will therefore 
be discarded here 9 . The configuration (B.l) with A" real p a is not the most general 
configuration in this case. An overall phase for all ip a cannot be removed anymore by 
t/yi (1) transformations. In contrast to the invariants p and Tj the properties of £ depend 
crucially on A^: £ is of order tp N . In addition there is an important difference between A" 
even and odd. For A" even the discrete symmetry ip — > — <p is part of SUl(N) x SUr(N) 
and is therefore respected by £. For A" odd £ is not invariant under this discrete symmetry. 
We also note that all invariants observe the discrete symmetry <p — > <p T and therefore also 

We restrict the discussion here to scalar field configurations with real diagonal eigen- 
values m a which are sufficient for deriving the flow equations for all derivatives of with 
respect to p and the T; (but not £). The discrete symmetry <p — > <p* of the potential 
implies invariance under <p/ j, — > — <p/ &. For real m a we can choose without further loss 
of generality the configuration (B.l) which respects ipi — > — </?/■ As an immediate con- 
sequence the mass matrix does not mix and <pj and (B.6) remains true. The mass 
matrices (B.4) and (B.5) acquire additional contributions: 

= V2— — (p a d a6T h <p c o 



9 We note that the sum uo 2 + £ 2 is proportional to det ip det <p^ and can therefore be expressed in terms 
of the invariants p, n. 
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+ ^ d 2 U k ( 5i 5r t 5£ 5r t 



didn \5p Rcd 5p a R b Sip a ^5cp Rcd/ 

(Plhj^J^ 8U k 5 2 j 
(dO 2 S<p% b Sp Rcd di 5ppp Rcd 
5 2 U k dU k 5 2 i 



+ 



8<f?5<p IcA " di 5pf5p Icd ' 
Here we have used S^/Sipi = for the configuration (B.l). Writing 

e = -L a — e b - b » ( Vaibl . . . VaNbN + < 6i . . . v * aNbN 



one has 









+ 


S( PRcd 




\°<Pcd 






(K 








\°<Pcd 


Wed 



with 



_ \ \ " ca 2 ...a N da 2 ...a N 

°<Pcd K N ~ !) ! a 2 ...a N 



^ \ ca 2 ...a N da 2 ...a N * 



oat 



a 2 ...apj 

for a diagonal configuration (p. Taking into account the reality of p> in (B.l) 
5£/5(pi = and 

= V25 cd n . 

For the particular configuration where all ip a equal p/N one finds 

JV-l 



and similarly 



^ 1 v 

5ppp$ ~ 5pf5pf (N - 2)! ^ e -3...e,^ e3 ... e ^ • • • V 



or 



5 2 Z 
5ppp<% 



(1 - 5ac) II Ve t 
eiy^a,c 

-(! - <w n ^ 
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In the following we will specialize our discussion to the particular configuration 



10 



( P<*=\jf) S <* (B.31) 

We restrict the discussion here to a potential linear in £ such that dU' k / <9£ = d 2 Uk/ d^dr 2 = 
d 2 Uk/d^ 2 = 0. The terms ~ S^/Sipji vanish in this case and the second functional deriva- 
tives read 



5 2 Z 5 2 £ / p 



N-2 
2 



= tf (! - Sac) (B-32) 



N-2 
2 



and for a ^ b 

5<ff*5<p1t~ 5ipf5ifff~ \n) ° ad ° bc - [ 66 > 

The mass matrices for the diagonal and off-diagonal ip a b as well as for ip^ and </?/ 
remain decoupled. For the diagonal part one finds 

ML = U' k 8 ac (B.34) 



jY 2 



The matrix Mj has N — 1 eigenvalues 



« = f (B.36) 



and one eigenvalue 



N-2 



(M?) 2 = U>-(N-1) (£) 3 ^ (B.37) 
whereas M| has TV — 1 eigenvalues 



N-2 



( M o\ 2 _ Tr/ i 4 P ^ (p\-^dU k 



and one eigenvalue 



AT-2 



(M|) Z = ^ + 2^p+(iV-l)(^) 2 ^. (B.39) 



10 



The reader should not get confused by our use of the symbol p for two different purposes — once 
for the invariant p = t\ and also for the field configuration p = Na 2 . Quantities like U' k always denote 
derivatives with respect to n at fixed T2 and £. On the other hand, for the configuration (B.31) the 
invariant £ becomes a nonvanishing function of a 2 — p/N . 
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The mass eigenvalues for the off-diagonal fields are 



K.) 2 = « + (^rf ■ (B.43) 



JV-2 



At the potential minimum one has U' k + (&f) 2 ^ = and we observe iV 2 — f Goldstone 
bosons ("pions") (B.36), (B.41), (B.43), one less than for dUk/d^ = 0. The scalar which 
has acquired a mass due to the ^Ta(I) violating term in the potential — the "r/-meson" 
— has positive mass squared (with dU^/d^ < 0), (B.37): 



There are also N 2 — f scalars in the adjoint representation of the unbroken diagonal 
SU(N) which have mass squared (M£) 2 , (B.38), (B.40), (B.42): 

r. JT N — 2 

(M° R f = Ml + lM| , Ml = ^^(Po.r, = 0, f = 2 (£) " ) . (B.45) 

Finally, there is the singlet ("radial mode" or V-field") with mass given by (B.39). 

We observe for p = that M 2 vanishes and M| vanishes for N > 2. In the symmetric 
regime one therefore has 2iV 2 real scalar fields with mass squared U k (0). This is differ- 
ent for N = 2 where M|(0) = —2^(0). In the symmetric regime the complex (2,2) 
representation decays in this particular case into two irreducible (real) representations 
with mass squared U' k (0) + ^(0) and U' k (0) — ^(0). In an obvious notation with Pauli 
matrices r k we can write the complex (2, 2) representation as 

V? = \ (a - irf) + \{a k + ™ k ) r k ■ (B.46) 

The fields (a, 7i k ) form an irreducible representation — they represent the standard linear 
(T-model [4] with scalars in a real four-component vector representation. In the standard 
linear cr-model the fields (?/, a k ) are omitted 11 and the scalar potential can be expressed 

11 This is a self-consistent truncation for N — 2, since (r/' , a k ) form a separate irreducible representation. 
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with the help of only one invariant p = |(cr 2 + n k 7ik). In terms of (cr, n k ) the Yukawa 
coupling (A. 6) with imaginary h reads 

^Yuk = ih a 1p (cr + i-K k T k m 1p 

_ _ _ _* (B.47) 
h = 2ih a , h a = h a 

with 7 being the Euclidean analog of 7 5 . 

If we decide to include the fields (r)', a k ) we obtain for the invariants p, r 2 and £ 

p = \ (vr a 7r a + a a a a ) 

r 2 = (a a a a ) (n b n b ) - (a a n a ) (a h n b ) (B.48) 

£ = \ (7r°7r a - a a a a ) . 

Here we have defined four-vectors ix a = (71^,0") and a a = (a^rf). We observe that all 
invariants in (B.48) are also invariant under the discrete symmetry a — > —a which is 
equivalent to ip — > The invariant violating this symmetry is u = i(det (p — det ~ 
a a 7i a . The quartic invariant r 2 can be constructed from p, £ and u;. 



C Evolution equation for A2 

In this appendix we discuss in more detail the f3 function for A 2 in the SSB regime. The 
evolution equation (3.22) can formally be written as 



d_ 

dt 



A 2 



= ~\l 



d d q d J N 2 A 2 



+ 



+ 



(2ir) d dt [ 4 (Z^P) 2 
9(N 2 - 4) A 2 N 2 X 2 



(Z V P + A 2Po ) 2 2p 



1 



6A 2 (iA 2 + A!) 



Z V P + A 2 p Z V P 
1 



Po(|A 2 -Ai) L^P + 2Aip Z tp P + X 2 po_ 



(C.l) 



where acts only on the infrared cutoff in P, i.e. t| = jr^§t~§p- Alternatively, we may 
use the form 





1 f d d q 'd J N 2 A2 9(A^ 2 - 4) 



dt 



X, = 



2 J (2n) d dt 4 (Z v Py 



+ 



[Z V P + 



N 2 



A 



+ 



3A 2 + 12AiA 2 



(C.2) 



2 (Z V P)(Z V P + \ 2Po ) (Z ip P + 2X 1 p )(Z^P + X 2 p ) 



AA 



which is close to the t-derivative of the contribution from the perturbative one-loop 
graphs, with p ~dependent masses and the infrared cutoff in the propagator taken into 
account. For d = 4, p — > 0, P = q 2 and omitting the t-derivatives one recognizes the 
one-loop correction for A2 [16, 17]. 

One may ask to what extent the (3 function for A2 depends on the choice of the 
configuration (3.13). For this purpose we have also evaluated dUk/dt for a configuration 

2 Po . 2 Po 2 Po r \ m 0\ 

¥1 = + e , V^2 = ^ - e , ¥a = for a > 3 (C.3) 

with infinitesimal e and 

From the contribution proportional to e 2 in dUk/dt one obtains the flow equation 

iL\ -I / A ^ / N(2N - 3) A 2 (2iV 2 + 5iV - 18) A 2 

^ 2 "2i(27p5t| 2 (Z^P) 2 + 2 (Z V P + A 2 p ) 2 

iVAa 3A2 + 12AiA 2 



+ 



(Z^P + A 2 p ) (Z V P + 2A lPo ) (Z^P + A 2 p ) 

(C.5) 

We observe that for iV = 2 the configuration (C.3) is equivalent to (3.13) and the evolution 
equations for \ 2 and X 2 , (C.2) and (C.5), respectively, agree. For N > 3, however, the two 
equations are not identical and this is the reason why we have chosen the symbol \ 2 in 
(C.5). The difference is related to the truncated higher invariants Tj, % > 3. The validity 
of the evolution equation (C.2) corresponds to a definition of the higher invariants where 
all Tj vanish for % > 3 for the configuration (3.13). Inserting (3.13) into (A. 13) we find 
T3 = — ( ' N ~ 1 ^~ 2 ^ (r 2 ) 2 and therefore (B.3) r 3 = 0. On the other hand, we find for the 
configuration (C.3) that f 3 vanishes (but not f 4 ) which leads to 

r 3 = (r 2 ) 3/2 . 

A truncation is well defined only once we specify the exact definition of the higher in- 
variants which we omit. This gives then a unique definition of the running of coupling 
constants. For N > 3 the difference in the results for k — > between the use of (C.2) or 
(C.5) can be taken as a rough estimate of the uncertainty due to the truncation of the 
higher invariants. 



D The infrared cutoff for fermions 

The choice of the infrared cutoff for fermions is not completely obvious for various reasons. 
First of all, chiral fermions do not allow a mass term. Since we want to remain consistent 
with chiral symmetries (a necessity for neutrinos, for example), the infrared cutoff must 
have the same Lorentz structure as the kinetic term, i.e. RkF ~ l^Qn [22]. On the 
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other hand, for q 2 — > the infrared cutoff should be ~ k, e.g. R kF ~ ^ / v^- The 
nonanalyticity of a/? m ay then be a cause of problems. We will develop in this appendix 
a few criteria for a reasonable infrared cutoff and finally propose one which seems suitable 
for practical calculations. 

First of all, the fermionic infrared cutoff term A k S F should be quadratic in the fermion 
fields. We next require that A k S F should respect all symmetries of the kinetic term for free 
fermions. This includes chiral symmetries and Lorentz invariance. (Gauge symmetries 
may be implemented by covariant derivatives in a background gauge field [28] but this 
if of no concern in the present paper.) The symmetry requirement implies in a Fourier 
representation 



A k S F = Q^TfMZ^k) (if a r F (^) V a (q) = ttqZ^r^ (D.l) 



where 77, 77 are Weyl spinors, a, a denote spinor indices, fa = q^ and cx M is the restriction 
of 7^ for left-handed Weyl spinors (with a suitable choice for right-handed Weyl spinors). 
We also have omitted possible internal indices labeling different Weyl spinors (for conven- 
tions see [22]). The wave function renormalization is chosen for convenience such that 
it matches with an approximation for the fermion kinetic term QjjZ^r] in T k . (Typically 
Z^p is diagonal in the internal indices but different fermion species may have different 
wave function renormalization constants.) The third condition requires that A k S F acts 

2 

effectively as an infrared cutoff. This means that for k — > 00 the combination Z^r F {%) 
should diverge for all values of q 2 . This divergence should also occur for finite k and 
q 2 /k 2 — > and be at least as strong as {k 2 /q 2 ) 1 ^ 2 . As a fourth point we remark that T k 
becomes the effective action in the limit k — > only if lim A k S F = 0. This should hold 

k— >0 

for all Fourier modes separately, i.e. for h.mZ^,(k) = const, one requires 

Sa" 0?) * = ■ (D2) 

Even though on a torus with longest circumference L the minimum value g^ in = n 2 /L 2 
does not vanish, we want a smooth limit to infinite volume and request that (D.2) also 
holds in the limit q 2 — > 0. Together with the third condition this implies exactly 



q 2 /k 2 ^0 k 2 ' 

The requirement (D.3) implies a smooth behavior of R kF for q 2 /k 2 — > 0. On the 
other hand, the nonanalyticity of r F at q 2 = may be a source of worry for practical 
computations. A careful choice of r F is necessary in order to circumvent this problem. 
First we note that r F appears in connection with the fermion propagator from Y k and we 
combine 

Z^+Z^rp = Z^F 

(D.4) 

P F = q>F 2 = q 2 {l+r F f . 
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Up to the wave function renormalization, Pp corresponds to the squared inverse propaga- 
tor of a free massless fermion in the presence of the infrared cutoff. We will require that 
Pp and therefore F 2 is analytic in q 2 for all q 2 > 0. A reasonable choice which we will 
employ in the present paper is 



Pf 



l-ex P {-£}' 



(D.5) 



E Scalar wave function renormalization 

In this appendix we provide some details for the calculation of the flow equation of the 
scalar wave function renormalization Z^^{k). We start with the scalar contribution to 
(4.5). It is convenient to decompose the complex fields ^> ab and (p ab into their real and 
imaginary parts: 



Vab(q) = -7= [<PRab(q) + ^iab(q)] , <plb(q) = ~m \vRba(q) - i<piba(q)} 



V2 



V2 



(E.l) 



For the configuration (4.1) and the ansatz (3.3) it is then easy to see that T Sk is block- 

(2) 

diagonal in the indices R and /, i.e. it decomposes into a block which contains 

(2) 

only functional derivatives with respect to the (pn a b and an analog block T Ik . We may 
therefore use the expansion (4.8) separately for T^ k and T$ . The first term on the right 
hand side of (4.8) is independent of 5ip, Sip* and thus can not contribute to the anomalous 
dimension. The part of the second term proportional to SipSip* is independent of Q and 
hence does not contribute either to the evolution of Z 9 ^(Q). Hence, we are left with the 
Q-dependent part of the third term whose contribution proportional to 5ld5lp* is obtained 
by keeping all terms of Ar}^ /A , which are linear in A(q, Q) = 5(p5(q, Q) + Sip*5(q, —Q): 



ay 



(2) 
Bk 



= A(g/ q ' Q) {(B + D)[6 ad E d , + 6 d £ ad ] 

ab,cd ip 



AT 



ab,cd 



(q, q') 



+ ^(C-B-D)[6 ab Z cd + 5 cd Z ab \} + 0(A 2 ) 
&(q'-q,Q) 



{(B + D) [5 ad E cb + 5 cb E ad ] 



D[S ab i: cd + 8 cd E ab }} + 0(A 2 ) 



(E.2) 



with 



From 



r 



(2) 
Rkfi 



ab,cd 



B = \m 2 , C = pU'^ D = ±M 2 



(q,q') = {(Z^ k (q)q 2 + U' k ) SaJbd + B (5 ac 5 bd + S ad S bc ) 



;e.s) 
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,(2) 
lk,0 



ab,cd 



+ ^(C-B)5 ab 5 cd }5(q, q ') 
(q,q') = {(Z^q)q 2 + U k )5 ac 5 bd + B(5 ac 5 bd -5 ad 5 bc )}5(q,q') (E.4) 



we obtain 



B 



/ F (2) D \-i , ,v / A(g) + £ ... 



2 C-5 
TV [/1(g) + 2C] 



[ r ^o + R k ) ab jQ,<l) = [ A{q) + 2B 5 ac 5 bd + A{g) - ^ a Ae f (E.5) 



where 

= ^(g)P(g) + ^ . (E.6) 

Now the traces can be evaluated and one finds for the scalar contribution to the scalar 
wave function renormalization 



with 

fk(Q) 



d 

-Q t Z v , k {Q) 



1 r d d q 
2~p& J (2^p 



= -%KM [fk(Q) - (°) 

s 



;e.t) 



d 



2pUj 



dt [Z Vtk (q)P(q) + U' k - j<Ml\[Z^ k {q + Q)P(q + Q) + U' k + 2pU>> - 



jn N-1 M 2 



and 



N 2 -4 



M 2 T + f Mf 



2 [Z v , fc ( g )P(g) + - iM|][^ ifc (g + Q)P(q + Q) + U' k + M* + ^Mf] 



Ml 



^M 2 



[^,fc(?)^(?) + ^ + P^(<Z + Q)P(q + Q) + U' k + Ml + ±Mf 



Ml = 



4p <9£4 
N -ldr 2 



Mf = -N yy 



p\ E 2 A dU k 



(E.8) 



(E.9) 



0£ 



Here we used the fact that P(q) and Z^^q) are actually functions of the invariant q 2 . 

We turn next to the fermionic contribution to the scalar anomalous dimension. Anal- 
ogously to the scalar contribution we may split Fp\ into Tp k0 and Ar^ which contains 
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all dependence on the scalar background fluctuation Sip. For the truncation (3.3) of the 
average action we obtain 



(41,0 + ^)^ M') = 



Zh(q)P F (q) + K(q)<p 

(Arp k ) ab (q,q') = h k (-q,q')&(q'-q,Q)^ 



$abS(q,q') 



(E.10) 



where we used the abbreviation h k (q) = hk(—q, q). Using an expansion similar to (4.8) the 
fermionic trace can be evaluated and we obtain its contribution to the scalar anomalous 
dimension by collecting all Q-dependent terms of order 5(f5(p*. One finds 



with 

fk(Q) 



d 



^(<7 2 + <?Q) h 2 k (-q,q + Q) 

Z^ k (q)[l + r F (q)] Z^ k (q + Q) [1 + r F (q + Q)} 

[Zl k {q)P F {q) + Ul(q)p] [ZlJq + Q)P F (q + Q) + ±)? k (q + Q)p] 



(E.ll) 



X 



d_ 

dt 



(E.12) 

and we made use of the identity h k (—q,q') = h k (—q',q). 

The fermionic contribution to i]^ obtains from the term linear in Q 2 in f Fk in close 
analogy to [23], and similarly for the scalar contribution, where we recover the results 
of [10] for A 2 = v = 0. Taking p at the minimum of the potential and neglecting all 
momentum dependence of the Yukawa coupling and wave function renormalizations we 
arrive at 



4 — -K 

a 



2Ai - 



N — 2 v ( k 



i 2 



2 N \N 



m? 2 (0,2«Ai —v 



N-2 
2 



;%) 



+ 



+ 



a 2 - 



N -2v ( k 



i 2 



N 2 -A 
2 



N \N 



V { K 



m 22 (K\ 2 + u[ — 



K 



N ( K 



i 2 



A 2 + - - 



N \N 



m£ 2 (0, k\ 2 + v ( — 



K 



;%) 



+ 2^N c h 2 mf )d {^h 2 - m ) . 

(E.13) 

The threshold functions and mf >d are defined in (4.11) and (4.12), respectively. 



49 



F Yukawa coupling and fermion wave function renor- 
malization 



We will give here details on the derivation of the evolution equations for the fermionic 
wave function renormalization constant and the Yukawa coupling as defined in section 5. 
We use here the truncation (3.3) but neglect for simplicity the dependence of the Yukawa 
coupling on the scalar momentum and also the momentum dependence of Z v . For the 
scalar fields it proves useful to introduce the following linear combinations: 



Via = y= (Vaa ~ V*aa) 



(F.l) 



and for a ^ b 



VR+ab = 


^ (Vab + Vba + V*ab + Vba) 


VR-ab = 


^ {Vab ~ Vba + V*ab ~ Via) 


Vl+ab = 


% 

~2 + ^« ~ ^ab - Via) 


Vl~ab = 


% 

— {Vab - Vba - V*ab + Via) 



(F.2) 



Using collective indices i,j G {Ra, la, R + ab, R — ab, I + ab, I — ab} this yields 

5 2 T k 



5 V *(q)5cpi(q') 
with scalar mass matrix 



(Z^b) + (M 2 )}) (2n) d S(q - q') (F.3) 



F) 2 TT 

(M 2 ) i = k . (FA) 

This matrix is discussed in appendix B and does not mix the various fields (F.l), (F.2). 
The mixed functional derivatives of Y k with respect to one scalar and one fermion field 
are easily seen to be of first order in the fermion fields. Since we are interested in the 
term bilinear in vp, we can take to be infinitesimally small and split 

if = rf ) + Aif (F.5) 

(2) (2) 

in such a way that all -^-dependence is entailed in AT k . Hence T k is given by (5.3), 

(F.3) and AT k 2 ^ is determined by the mixed scalar-fermionic functional derivatives of r*.. 
Using an expansion similar to (4.8) we arrive at the following evolution equation for the 
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bilinear fermionic part T^l of the effective average action: 



dt k > 2 ~ 



J_ rd*Q_ r d d g T 2 g + Q d_ f ^(g)t(Q)^l+r F (g)]f(Q) 
2N J (2k)' J (2k)' * 2 ^1 Z$J q )P F (q) +T? k (q)<? 



(27v) d J (2tt) 
N 2 -l 



X 



+ 



+ 



X 



+ 



N 2 ~ 1 



Z v P(q -Q) + U' k - ±M* Z^P(q - Q) + U' k + M* + ^Mf 
1 1 



+ 



Z v P(q -Q) + U' k + 2pU' k > - ^M 2 Z^P(q - Q) + U' k 
h k (q)^ a (Q)lipa(Q) 



N-l M 2 



N 2 



Zl k (q)P F (q) + hl(q)<? 

{ N 2 -l 

\Z v P{q -Q) + U' k - JM 2 ~ Z v P(q -Q) + U' k + M? + ^Mf 

1 1 

+ 



Z v ,P(q-Q) + U k + 2pU>>-^M 2 • Zip P(q-Q) + U' k + ^M 2 )l 

(F.6) 

with M 2 = and Mf = -N (f) 2 It is easy to extract from (F.6) the 

evolution equations for the Yukawa coupling and the fermionic anomalous dimension. For 
the Yukawa coupling, renormalized according to (5.2), one finds for p at the minimum of 
the potential 



[d - 4 + 27fy )fc (Q) + 7fo fc ] h 2 k (Q) 



+ 



x 



+ 



k 



4-d 



N 



MQ) 



d 



(2nY hk{q)hk{ 2 } Z,, k (q)Z,, k (Q) Z * 
N 2 -l 



dt\P F (q) + ^hl(q)K 
1 



Z v P(q -Q) + U' k - jfMl 
N 2 -l 



(F.7) 



Z v P(q -Q) + U' k + Z v P(q - Q) + U' k + M* + ±M 2 

1 



Z„P(q ~Q) + U' k + 2pU' k > 



II N-l 
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Analogously we obtain for the fermion anomalous dimension 



k'~ d f d d q qQ h2( q + Q, Z\ k (^)Z v fl ( [l + r F (q)} 



Vl>,k(Q) = o~aF / 7^ h , 



2N J {2txYQ 2 ky 2 1 Z^ k (q)Z^ k (Q) dt \ P F (q) + f h 2 (q)K 



( N 2 - 1 1 



+ 



+ 



V(« -Q) + U' k - ±MI Z^P(q -Q) + U' k + ^M 2 
N 2 -l 



Z^P(q -Q) + U' k + M 2 + ±M* 
1 



Z v P(q -Q) + U' k + 2pU» - ^M 2 ) I ■ 

(F.8) 

In the following we do not discuss further the momentum dependence in the fermionic 
sector. We are interested in the evolution of rj^ = %,fc(0) and h 2 = h\(Q). Neglecting 
the momentum dependence of an d hk on the right hand side of (F.7) and (F.8) and 
taking the limit Q 2 — > we obtain 



h 2 = [d-4 + 2 % + ^]/i 2 --^ 4 {(A^ 

1 N N ~ 2 

+ l{ F i B)d (^Kh 2 ,e+ —v (^j 2 ; 77^, 77^) 



(iv 2 -i)C )d (^ 2 ^ + ^ 2 + ^(^) 2 
^ = ^^ 2 {( JV ' 2-1 ) m S B)d (^ : ' i2/c ' e ;^'^) 

JV-2 

+ ( J V 2 -l)m< F 2 B) ' i (-fc 2 K ,e + KA 2 + i'(^) 2 ;%,'!».) 

The threshold functions and m^^if are defined in (5.5) and (5.7), respectively. In 

the symmetric regime one has to set k = in both expressions, whereas in the SSB regime 
e = 0. 



(F.9) 



(F.10) 
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G Threshold integrals 

In this appendix we provide explicit expressions for the threshold integrals introduced 
in the previous sections which are suitable for a direct numerical integration. For this 
purpose we introduce the dimensionless integration variable y = q 2 /k 2 = x/k 2 . Using the 
abbreviations 

P = p{v) = k~ 2 P(x) , p = p(y) = , etc. 

p F ee PF (y) = k~ 2 P F (x) , p F =p F (y) = ^01 , etc. (G.l) 



we find 



_ / x Pf{v) n • / \ dr F (y) 

r F = r F {y) = y — 1 , r F = r F {y) = — , etc. 

l*( w ) = n r dyyi~ \ P ~ yp (G.2) 

ny 1 Jo [p + w] n+1 v ; 

ii(w) = |/; d „i- lF ^_ (G . 3) 

■ «"> 

l™{v,) = U r d yyi I^L+A (G.5) 

l i una ( w uW2) = J q dyy2- 1 (p-yp) 

x It r^r ^ + ? ml (G.6) 

{ [p + w 1 ] n i+ 1 [p + w 2 ] n2 [p + wi] ni [p + w 2 ] n2+1 J 



I r-oo d 



x t r^n ^ + t m ( g - 7 ) 

/■oo d I 

G>^ 2 ) = j o dyy^ [pF + Wi]ni[p + w2]n2 
x U[p F -yp F ] + n 2 [p-ypn 

I [PF + ^l] [P + W 2 ] J 

/f2 d K, W2 ) = ^.Tdyyi- 1 -. ^ t-t (G.9) 

ni ' n2 v ' ; 2 Jo yy [p F + w 1 }^[p + w 2 } n2+1 V ; 

= ^Tdyyi-. ^ferr 1] — r- (G.10) 
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x j n^jp-yp] | n 2 p\p-yp] | ] 



i f 00 

m d nin2 (w 1: w 2 ) = -y o rfyy 



P 



x 



[p + w l ] n ^[p + w 2 } n2 
1 [p + ^i] [p + w 2 ] 



[pi^ + w] 3 

x S [?V + yr F ] [p F - w] - yf F p F 



p[l + r F ] . ' 

- r F 



p F + w 

p F — 3w 



m[ F)d (w) = I dyyi +1 - — T 

J \p F + w\ 



p F + w 

p[l + r F ] 
p F + w 



-r F 



j . r 1 . P F ~3W 

x <r F [p F - w\ - r F p F ■ 

I p F + w 



|[i , 1 I P\p-yP] . ~\ . . ( Znxpp 
x <[l + r F \\n 2 ■ \-yp]-ypr F 



p + w 2 J XPf + W! 



2 Jo [p F + wi] ni [p + w 2 ]" 2 

p[p-y] , l 



x < n 2 h 1 - p 

I p + U>2 



,(FB)d 



d 



r F p / 2riip F 
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